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9m  prtMSt  rapoort  It  a  pwt  at  the  r*Mareh  progxHi 
OB  bouadazy  liyar  prbblMM  purauad  at  tbe  draduata  Ckdiool  of 

Aarooautloi^  BigiiMazlag  «««<*»»  tha  provlalona  of  eoatxaot  AF  49(638)» 
9)is  eoBtraet  la  flnaaeed  Tigr  tba  italted  Stataa  Air  Force  Office 
of  Sdeatlfle  Biaaazch  and  la  mooitorad  by  the  Neehaales  OLvialon. 

9m  isteat  of  thla  lanreatlflatloa  baa  been  to  eaat  on  aoaw 
of  bouadazy  layer  flow  at  atagaatloa  polnta  In  flow  of  air 
at  apeeda« 

9ila  lanreatlgatioii  imm  propoaed  by  fttxfeaaor  Hteholaa 
B9tt  aad  the  lareatlgator,  Mr.  Leaard,  haa  woriced  under  lha  direetloii 
of  nrofeaaor  Bott  and  nrofeaaor  S.  E.  Laa.  9ie  research  carried  out 
under  thla  oontraet  la  under  the  general  dlsreetlon  of  Rrofeasor 
V.  B.  8eara«  Director  of  the  School. 


ERRATA. 


P.  16,  Equation  (l.23)i 

Aid  factor  of  (l4>r))  on  rlght»hand  side  of  equation,  vis.i 

P.  21,  Equation  (2.8)t 

6th  line;  omit  R  in  ienominatepi.  of  1st  ani  2nd  terms;  change  R 
in  denominatbift:-  of  3rd  and  ^th  terms: 


4.  5«  4-  _ hu>| 

7th  line;  omit  R*-  in  denominatiMJn  of  1st,  2nd,  l|th  ani  5th  terms, 
viz.: 


I  i 


8th  line;  change  density  derivatiTe  in  last  term: 


|v-j 

H  l.» 


+  V«/  J 


♦y.,A  J'+Yv^ 


P.  22,  Equation  (2.13): 

2ni  line;  change  2nd  term  from  X  to  ^  derivative: 


P.  65,  Uth  line  from  bottom: 

Insert  j  ,  viz.: 

.  terms  with  the  J  subscript 


P.  65,  2nd  line  from  bottom: 

Insert  K,  and  ,  viz.: 

.  of  proportionality  and 


P.  9^,  Equation  (F.ll): 

In  right-hand  aide  of  2nl  and  3rd  expressiona,  replace  h  (enthalpy) 
by  jju  (viecoaity),  viz.t 

&.«»  •  ^  Jw  (■» 

P.  95,  Equation  (F.12)» 

Laat  line  refers  to  function  ,  viz.: 

jr' 

('ll  — ►  (»v“  I  )‘^ 

P.  103,  Equation  (H.2): 

Ra.  on  right-hand  aide  has  subscript  I  ,  viz.: 

. 

A“ 

P.  105.  In  line  above  equation  (H.?),  reference  is  made  to  (H.3)  instead 
of  (h.3)f  viz.: 

(H.3),  are  presented  . 

P.  1C6.  In  definition  of  ,  replace  V)  and  V  by  yiu  and  0  ,  viz.: 

^  shear  parallel  to  surface  ^  *  M. 


ABSiaUTP 


StMdy^  vlsoou*>  tw»>dlMraaioBAl  nd  aartalljr  qnMtzie 

■tacMtloi^polnt  flows  of  a  gu  art  ooBsidtrsd  for  tte  aasa  idam  tbs 

atyBOlds  Buabar  is  too  low  for  tbs  arolloabllity  of  tbs  elasslMl 

boandaxyblsyar  tbsoiy.  It  Is  assuasd  that  Iba  lflw»4sBSltgr  gas  Is  still 

a  eontlouoas  fluid,  pezslttlag  tbs  uss  of  tbs  ftnrlax^Btdces  and 

assodatsd  sgaatlODS  as  tbe  basis  of  tbs  proUaa.  flw  offsets  of 

low  fteynolds  ausiber  are  dstemliksd  bgr  ^nOjrlag  sa  sbjsmIob  proeedure 
.  (20) 

(slaUsr  to  Lsgerstroa  aad  Cola's^  '  la  taasa  of  a  panMStar, 
esssotlally  ,  to  tbs  flui*>d]nHaisal  •fMtlons..  OIm 

bl|di«*'^<nder  eqtiatloas  la  this  saqpansloa  axe  tbs  boualaxsMagrsr 
equations;  tbs  osxt»ordar  equations,  idiiah  tbexefors  Ixrol've  flrst>order 
low>BS3mal4»>xiuBiber  coxTsetlon  tszas  to  tbs.bowaaaxgr-laver  qjoaatltles, 
are  presented  sad  discussed  la  detail,  togetber  wltb  tbe  appropriate 
boundary  oonditlons.  Ibe  boundary  eondltlons  axe  of  tso  tgfpes;  at 
tbe  sail  tbey  are  derived  fron  tbe  klaetie  ttasosy  of  gMes,  fbr  fron 
tbe  sail  tbe  flow  oust  "aerge"  Into  tbe  Ixarlsold  s«drrtioii. 

It  Is  shown  that  tbe  foUowlig;  qaaatitleB  axe  aeeessaxy 
to  define  tbe  Isvlseld  flow  near  tbs  stagnation  point:  tbe  ataguatloa 
properties  of  tbe  gas,  a  veloeity  gradient,  a  nose  xndltts,  aad  a 
Tortlelty  paxeaater.  3be  latter  Is  present  In  tbe  aadally  sjnnHtxle 
ease  oinly;  It  defines  tbe  slope  of  tbs  Inrlseid  sbsar  flow  near  tbe 
stagnation  point  of  tbe  axially  a|nBetxf.e  body.  It  is  dunn  that 
furtber  generalisation  of  tbs  Ixarlseil  flow  by  eonaidexlai  additlanal 
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p«z«wt«n  MeasMkxy  to  dtflao  tha  flaw  otrar  iBqpar  raglou  ia  tte 
vleiaitgr  of  th«  stacnatioa  polat  la  aot  BaeaBaax7f  taoanaa  aagr 
additional  paroMtara  vUl  aot  affaot  tba  llrBt>OKdar  eoanaetlona 
oonaldarad  la  tba  aaaljraia  of  tba  rLaaoua  Haw. 

Oia  raault  of  tba  aaalyala  la  tbat  tba  fdUoalaK  eorxaotloa 
affaata  to  bouadaz3P«layar  tbaocgr  ara  daflnad:  euxvatara  affaet« 
dlaplaeaaant  effect,  veloeltgr  allp  and  teagpexatura-Jiav  effect, 
Tortleltgr  effect  (thla  laat  one  for  tba  axlaUor  agraaetzle  oaac  only). 
9m  aagnltuda  of  tbaae  effeeta  dapeatbon  tba  foUoadag  pazaaartera, 
zeapeetlTeljr:  tba  ratio  of  tba  bouadaijr^lagrar  tblohaaaa  to  aoae 
radliia,  IIm  obaaga  In  atagnatio»>polat  Talocltgr  gzadlaat  ffoa  tba 
Isvtacld  value  due  to  the  dlaplaeeaeat  effect  of  tbe  bonadaxy  layer 
around  tbe  body,  the  ratio  of  the  aaan  free  path  to  the  boundaryw 
layer  tblekneaa,  and  finally  tbe  ratio  of  tbe  dope  of  the  larlaeld 
.8taear>flov  profile  (in  axially  lyaatrlc  flow  only)  to  tbe  acremge 
dope  of  tlM  boundary>layer  velocity  profile. 

Sone  appUcatlona  are  diacuBaed,  In  particular  tba  oaae  of  a 

a 

blunt  body  flying  through  the  ataoapbare  la  oonaldarad  In  detail. 
Beal^gas  propertlea  are  uaed  to  calculate  the  cxpanalon  paxeaMtera  for 
thla  caae;  tba  reaulta  are  plotted  in  chart  fom.  9m  region  of  beat 
appUoablU-y  of  tbe  expaadcn  procedure  la  vtaare  the  e^paadon 
paraaetera  are  leaa  than  1  and  ara  about  the  aana  order  of  ancattude. 
9>la  oeeura  In  tbe  fll^t-apeed  range  of  Ihcb  2  to  7« 

ICnexlod  adutlona  of  the  efuatlona  are  than  praaentad  ualng 
tba  propertlea  of  undlaaoelatad  air  oorraapondlag  to  thla  apaad  range 
(perfect  gaa,  conatant  Praadtl  anbbcr,  variation  of  dpadflc  beat, 
vlacoalty,  boat  eonduetivity  dtb  poaera  of  abaduta  tenparature). 


BBBulta  art  tabulated  for  five  values  of  vail- to» free- atrasa  taagperature 
ratio;  exaag^a  of  lofie-asynoT ae-iwiabar  veloeltgr  aad  teapezatura 
profiles  are  given.  •Ban  effect  of  taaqperatura  ratio  oa  vail  dhaar 
aad  vail  beat-traaafer  rates  are  shewn  in  both  table  aad  graph  focas. 

Hie  behavior  of  the  shear  aad  hea'b>tznasfer  eorreotiaas  due  to 
veloeltgr  slip  aad  teBqpezature  jump  is  eq^lally  sigalfieaat;  the 

» 

results  elearly  indicate  that,  in  spite  of  the  ouch  wHller  aeaa  free 
path  at  the  vail  for  stronglgr  cooled  boundary  layers,  the  rednetloa 
in  heat  transfer  due  to  this  effect  is  detezained  by  the  aeaa  free 
path  in  the  (hot)  Invlseld  flow  at  the  stagnation  point.  Oils 
conclusion  is  the  result  of  <n«i using  the  effects  of  variable  fluid 
properties  in  the  analysis. 

Qie  predicted  stagnatloiwpodLnt  heat-tzansfer  reductions  for 
a  cylinder  in  a  supersonic  alrstreaa  at  low  Beynolds  nuaibers  are 
coapared  to  the  experimental  results  of  lewfik  aad  Qledt^^^'^^^^. 

Both  theozy  aad  experiment  indicate  reductions  in  heat  transfer  at  low 
Beynolds  nuaibers;  but  the  measured  reductions,  vhile  agreeing  in 
trend,  are  considerably  larger  than  the  prediction*. 


iwaoDDOiiai 


Die  flow  In  the  vicinity  of  fonnrd  stegBaitloa  points  has  long 
been  of  qpeolal  Interest  to  aerodynsalclsts.  The  reasons  for  this 
interest  are  both  theoretical  snd  practlesl.  !Zlheoretlcally>  the  flew 
Impinging  on  sn  "Infinite"  flat  plate  provides  one  of  the  fev  exact 
solutions  of  the  ]favle]> Stokes  equations;  furthemore  the  solution  of 
boundary^lsyer  flew  around  bodies  alwys  "beglits”  at  the  stagnation 
point.  Practically^  heat^tzansfer  rates  are  usually  at  a  iwarlmnn  In 
the  stagnation  region;  furthemore,  stagnation  properties  of  the  flow 
are  often  the  easiest  to  neasure  by  aeans  of  various  probes  and  measuring 
Instruments.  In  addition.  In  recent  times,  the  reentry  proiblem  has 
generated  renewed  Interest  in  this  problem,  especially  in  the  low* 
Heynold8>l(umber  flow  regime.  The  purpose  of  the  present  Investigation 
Is  to  reexamine  this  problem  In  this  loM>>Itoynold»>nuid>er  flew  regime 
for  the  Bteady>flaw  case  and  for  the  plane  two>dimenalonal  and  axially 
syanetrlc  flow  patterns;  due  to  their  slmHaritles  a  unified  treatment 
of  the  two  types  of  flew  will  be  possible. 

loqillclt  In  the  aerodynamiolst  ’  s  solution  of  flow  problems 
Is  the  assumption  of  a  continuous  fluid.  Ihe  classical  eqiiatlons  of 
Navler^Stokes,  together  with  the  continuity  and  energy  eqjuations,  and 
additional  equations,  idilch  relate  the  themodynamio  and  transport 
properties  of  the  fluid,  fom  the  basis  of  continuum  aerodynamics. 

This  approach  will  be  maintained  In  the  present  Investigation;  though 
It  must  be  duly  noted  that  It  Imposes  a  very  definite  limitation  on  the 
applicability  of  the  results  on  gases  (which  are  principally  of  Interest) 
In  terms  of  a  Khudsen  nuoiber,  which  cannot  exceed  a  "reasonaUy  small" 
value.  What  this  limitation  means  can  be  Inferred  from  the  ULnetlo 
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theory  of  gases.  It  Is  veil  kBoim  that  the  properties  of  a  fluid  flow 
field  for  a  gas  eeu  be  obtained  If  the  (aoleeular)  reloelty  distribution 
function  is  given,  by  talcing  the  appropriate  statistleal  averages. 
Solutions  of  the  Mucvell-BoltBaBan  eqiiatlon,  Vhloh  is  the  conservation 
equation  for  this  distribution  function,  can  be  obtained  in  terns 
of  an  Iteration  procedure,  which  results  in  successive  sets  of 
“continuua"  differential  equations;  i.e.,  eqjuations  in  terns  of  the 
"locally  average"  quantities.  It  has  been  pointed  out  by  a  nunlber 
of  authros  (e.g.  Sheroan^^^^^  SohaSf  and  Cbaadtre^^^^,  etc.)  that  this 
iteration  procedure  is  roughly  equivalent  to  an  expansion  in  terns  of 
the  oean  free  path  between  aoleeules.  Ihe  successive  sets  of  equations 
are:  Euler's  (inviscld}  equations,  the  Havlerw Stokes  and  associated 
equations  (which  will  be  used  in  the  present  work),  the  Burnett  equations, 
and  equations  of  still  hifljuir  order.  It  can  be  seen  then  that  the  present 
analysis  ttegleots  the  Burnett  and  Idgher^order  terns,  vhidx  iogplies  that' 
only  "snail"  changes  in  local  flow  quantities  are  pemlsslble  over 
distances  equal  to  the  mean  free  path.  Actually  it  has  been  fmnd 
(Schaaf  and  Chandire  that  even  for  sqsm  specially  simple  and 

well  understood  prcblems  in  gasdynsnlos,  such  as  the  structure  of 
nonsal  shocks  and  the  propagation  of  hlgh»frequeney  sound  wanras,  where 
the  Burnett  texips  were  clearly  nonpnegUgible,  their  inclusion  in  the 
theory  gave  lees  satisfactory  agreaaent  with  experiaental  results 
than  the  theory  based  on  the  BavieivStokes  eq^ations  only.  Furthemore, 
the  Burnett  terns  Include  derivatives  of  higher  order  than  are  present 
in  the  navieivStokes  equations,  which  Indicates  the  necessity  for 
additional  boundary  conditions.  Ihere  is  no  agreenent  at  the  present 
as  to  what, if  eay,  these  additional  boundary  conditions  should  be.  A 
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more  detailed  dlscuBsloa  of  these  and  other  difficulties  ixiTolved  la 

a 

the  use  of  the  Buznett  equations  can  he  found  in  reference  ^6. 

Qie  present  woric  is  then  on  application  of  the  Savlez^Stokes 
equations  to  the  stagnation-point  prchlea  for  a  gas  vlth  knoim  thenm- 
dynanic  and  transport  properties.  A  perfect  gas  will  he  assumed^  with 
the  specific  heat;  vlscosity>  and  heat  conductivity  proportional  to 
arhitrazy  powers  of  the  absolute  temperature,  (ihese  assuBqptions 
regarding  the  properties  of  the  gas  will  he  Justified  in  another  section. ) 

Qie  classical  solution  of  this  problem  for  vezy  large  Reynolds  nunbers  has 

* 

two  phases:  first;  the  solution  of  the  Invlscid^flov  problem  with  slip 
around  a  given  body  is  obtained;  this  gives  the  location  of  the 
stagnatiaa  point  and  the  velocity  gradient  there.  !Ihe  second  phase 
completes  the  calculation  of  the  flow  by  applying  these  results  to  the 
well  known  solutions  of  the  stagnation-point  boundazy-layer  eqioations  in 
the  two-dimensional  and  axially  syasKtrlc  cases  respectively  (e.g.  Cohen 
and  Beshotko^^^,  Brown  and  DonoucdM^^^;  Howe  and  Mersman^^^^;  etc.).  It 
is  shown  by  Lagerstrom  and  Cole^^^  that  these  two  steps  can  be  looked 

I 

upon  as  the  first  steps  in  an  ejqpenslon  procedure  for  obtaining  approxi- 
matlons  to  the  solution  of  the  Hhvlez^ Stokes  equations  for  high  Reynolds 
numbers.  Uhls  procedure  consists  of  expanding  the  stream  function  in 
powers  of  a  Reynolds-number  psirameter;  (in  this  case  ^ 

in  tezms  of  two  parallel  series,  the  so-called  "Inneif  and  "outer" 
ejqpanslons.  Successive  tezms  in  the  two  series  are  solved  for  alteznating^ 
ly.  In  this  naaner  first  the  Invlscld  flow  around  the  body  1'*  obtained 
(the  first  tezm  in  the  "outer"  expaasioa),  then  the  coaq^ete  boundary- 
layer  problem  solved  (the  first  tezm  in  the  inner  expoasioa);  then 
a  correction  to  the  external-flow  foUcwi  (usually  due  only  to  the 


dloplaceneat  effect  of  tiie  bouadaxy  layer);  end  then  a  oorreetlon  to 
the  bouxtdaxy- layer  solutloa,  etc.  It  la  apparent  frcm  phyaleal 
conalderatlona  that  the  foxmal  aathematleal  procedure  developed  in 
thla  reference  la  rigorously  applicable  only  to  certain  very  apecial 
types  of  problems.  It  entirely  fails  to  account  for  such  \u)lversally 
present  phenoneca  as  turbulence;,  and  the  posslblllly  of  separation^  for 
example. 

In  the  foUawlng,  a  treataient  of  this  second  approadniation  to 
the  "Inner"  flow  for  these  two  t^pes  of  stagnatlon-polnt  flows  Is 
presented.  Cm  formal  mathematical  procedure  of  reference  20  will 
not  be  followed;  still  the  above  coanents  Indicate  that  a  complete 

f 

treatment  of  this  secoord  approximation^  which  is  essentially  an 
is^rovement  of  boundary- layer  theory  for  low  Reynolds  numbers,  would 
require  a  solution  of  the  second  approodmatlon  to  the  "cuter"  flow. 

This  Implies  that  In  addition  to  a  supeximposed  external  velocity 
gradient  of  arbitrary  ms^pltude  (which  Is  the  result  of  the  solution 
of  the  first  approximation  to  the  "outer"  flow,  l.e.,  the  Invlscld 
flow)  one  should  have  an  additional  arbitrary  external-flow  parameter, 
tdilch  influences  the  second  "lnner"-flow  appraxlmation,  and  which 
should  properly  be  the  result  of  calcxilstlng  the  Improvement  to  the  outer 
flow  due  to  the  displacement  effect  of  the  boundary  layer.  Obviously 
such  a  calculation  cannot  be  made  within  the  fremewoxlc  of  calculating 
the  flow  at  the  stagnation  point  only.  Cms  one  has  to  accept  the 
presence  of  an  addltiooiBl  arbitrary  parameter  in  the  problem,  due  to 
an  unleteimlned  "displacement"  effect. 


CHAPTER  I 


InvlBcld  Floy  near  the  Stagnation  Point 

Let  R  be  the  radius  of  cuz7*ture  of  the  body  near 
the  ata^natloa  point.  Then  the  "polar"  coordinate  system  of 
Appendix  A  (A. 6)  can  be  transformed  Into  the  conventional  boundary- 
layer  coordinate  system  by  the  trazxsformatlon: 

LJ.  r  -  R 

(1.1) 

\diere  defines  the  body  contour,  as  shown  In  the  sketch. 


Using  transformation  (1.1)  in  (A.8)  and  (A. 9)  the  imrlseid 
momentum  equations  can  now  be  written  down.  In  the  ^  direction; 


1 

1*1 


('*?)' 


i  +  n 

C'^  If 


(1.2) 
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aiKUiuQy,  in  the  x  dlreetlon 


f  >  'f'-  't''* 


5P.=o 


0-3) 


The  velocities  can  be  related  to  the  eoogpressible  stream 
function  by  tislag  (A.6)  and  (l«1)  in  (A.3): 


I 


u  — 


■r  =-. 


1-  . 


(lA) 


Since  flow  in  the  free  stxeam  la  ualfom>  and  the  flow  is 
steady,  the  "Iso-enereEici'ty"  condition  can  be  written  as: 


I 


'■'»  *-  iTTi^ 


vl** 


K  -  —  0 


(1.5) 


where  the  subscript  ^  refers  to  stagnation  eonditionB.  Tor  the 
assunptlon  of  teaperature-dependent  speeifle  heat  the  enthalpy  can  be 
expanded  about  the  stagnation  condition  in  a  Thgrlor  series,  as  follows: 


where  the  dots  slgnlf^friEKteilvatives  with  rsbpeet  to  tanpazntore. 


(1.6) 
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The  qmetzy  of  th*  prtHaajmh  and  tte  bounaugr  eondltlon  that 
vr  =.  O  on  th*  -body  «  O  eaa  b*  uMd  to  mite  down 

the  foUoMlng  eaqpanBlooa  for  the  various  qiaaatltiss: 

%-  Ss  ['  ^  a-olr  +-•>  +  +  ] 

T  -  +- a-Te^^^“*'+ 

P  '  P.  L'  <»•<  S,x^  Si-S— ] 

0.7) 

33m  perfect>-gas  lev  can  ho  used  to  ralato  the  prasaun, 
deaaltgr,  and  tonperature  erpaaslona: 


O-V,  -  ^'-Tl  ‘^r| 


a,p^  —  ^Ti  ®'ri  ■*■ 


Wf  *  Wt  Si  ^  S,  +  •-»»  V, 


n 


0.8) 


11 


expaaaloas  (l  .6)  and  (1 .7)  wlLL  b«  aUbBtltutad  Into 
•gMAtians  (l.2)>  (1.3}>  aad  (I.3).  Coeffielenta  of  like  powers  of  tbe 
Indapeadent  vmzlableB  will  then  Im  eq|Mkted>  gLvlag  a  suoeessloa  of 
additional  relations,  similar  to  (1>8).  Ibe  results  of  tMs  procedure 
will  be  that  many  of  the  coefficients  In  expansions  (1 .7)  will  be 
expressed  in  tens  of  a  smaller  mniber  of  independent  ones.  Oils  will 
Show  the  truly  significant  independent  parawters  neoessazy  to  describe 
the  Inelseld  stagnatlocwpolnt  flow  that  we  are  considering.  Ibe  details 
of  the  procedure  follow.  Substituting  the  e^gpenslons  (1.6)  and  (I.7) 
Into  the  energy  equation  (1.3)>  coefficients  of  the  tern  give: 

-  thus 

at,  -  0 

using  this  result,  the  ^  ,  and 

^Si  -  0 

-*■  a  5*  0 

.  .  S/-  A^c,  r  A,:  c,T,  ^  .  a  c, J,  -  0 

froB  which 


(1-9) 

terns  give: 


I 


(1.10) 
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OL 


Tfi. 


a  - 


-  -  ^  <r. 


(1.11) 


a  '  (1.12) 

A  b1  mi  1  ftr  substitution  Into  the  aoDentun  equuatlon 

(l  .2),  and  grouping  coefficients  of  the  constant,  ,  and  x** 

terns  gives: 

SiPsO-vv^  ® 

+  2.s,PjO.^^+  s,  p,  =.  0 


fron  uhich 


a. 


(1.13) 

■  ^Vs 

(1.14) 

.  5s 

9s  "R 

(t-15) 

finally  the  X  .  aonentun  equation  (1 .3}  is  used 
in  a  ejillar  vagr,  groupiag  coefficients  of  the  constant  and  tezau: 

'*■  5sPi  K  =  0 

t  -  2(k«)o-,  ^1+  2-  *  0 
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heace, 


(1.16) 


How  use  (1.8)  together  mth  (1.9),  (l.lO),  (1.11),  (1.13),  (l.ll»). 


sad  (1 .16)  to  get: 


a. 


r\ 


0 


(L 


- 


I 

=  a-cj,  <SL 

(1.17) 

Diese  results  aske  It  possible  to  detexniae  the  coefficleat 
of  the  second  tezm  in  the  ezqpeoisioQ  for  the  streu  fuactlon,  d.,  : 


Thus: 

for  m-O  cu,  =  -£^ 
for  n  =.  I  ou,  iB  axtdtrazy 

Aa  expression  for  the  vortieitjr  near  the  stagnaticii  point  can 
be  fouM  by  using  (1.1)  in  (A.  13): 


4- 


%  % 


\%r  "I 

1)  J 


(1.18) 
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Subatltutlag  uqpaaalon  (1.7)  into  (1.18)  tte  «qpc«sil(m  for  tb* 
vortieltgr  Is: 


A*  -  Aj.L2.o-.  +  J  ^  . 

Thus,  for  tha  two  Igrp**  of  sta0Mttioo>polBt  flow: 
for  n  «  O  il  *  CTcx.^) 

for  J1  - -A>.[iiou,  4Crt^)]a-A*[V'+‘j 

(1-19) 

wbera  the  shore  expression  defines  a  rortieUgr  puMseter  V 
Ooeffleleat  a. ,  can  then  he  written  as 


2-  (i.ao) 

This  result  can  now  he  used,  together  wltii  (1.8)>  (1.12), 
(1.1$)  and  (1.17)  to  detezaine 


1,  =  -  r„v  -  i  1 


(1.21) 


flie  shore  intezaedlate  results  san  he  siasHuri.xod  hgr  writing 


ejqpaasions  (1  .J),  using  the  dsrlrad  mloes  of  the  eoefflelaats: 
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4.=  '■'] 


(1.22) 


T.T.-£[(,..,^|^.|'-(i-v).S....J 


(1.24) 


(1*25) 


nroB  these  eaq^resslotis  the  pueaeters  for  the  qnestrloel 
Imrlseld  stagostlooppolnt  flctir  eaa  he  aseertalaed.  Veloeilgr  giedleat  A 
Is  a  fuadseental  pazeaeter;  It  is  dstezaiasd  hor  solvlag  the  InrLseld 
flow  arduad  the  entire  body.  It  will  be  of  the  order  of  aacpaitude 
of  ftn^stresB  Tsloolty,  U  ,  divided  bQr  a  "aicpdfleaBt  body  size," 

L  t  perpendicular  to  the  flow  direetion.  Tvo  sore  paxaoeters 
appear  in  the  streaa  function;  nose  radius  R  ,  and  vortielty 
parsawter  V  .  It  has  bean  pointed  out  prevloualy  (e.g.  Bott 
sad  Lsaard'  ')  that  there  is  an  lapoortaat  differeaee  in  this  respect 
between  tw>dlaBnBiQoal  and  axially’  qnaaetzle  stagaatloBppolnt  flows;  the 
latter  adaittlag  a  vortical  teza  but  showiag  no  effect  of  nose 
curvature  to  the  order  of  terns  that  axe  considered  here  (i.e.,  )* 

idMzeas  the  twodlaensiooBl  flew  is  irrotationsl,  but  has  a  curvature 
tern.  Stagnation  thenaodynealc  properties, 
and  au:e  the  additional  azUtrazy  pazeaeters.  Bud 
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lubMgueat  terms  la  the  eaqpaaslon  (!•«•/  bii^xworder  terms  la 

end  t^  Included  la  (1 .22)  through  (1 .8$)  )  heea  considered) 

addltlonsl  parsmsters  vould  have  appeared,  flhese  panmeters  will  he 

azhltniy)  or  IdentUtlahle  as  the  rate  of  ehsage  of  the  nose  mdlus 

vlth  )  etc.  The  possible  effect  such  additional  arhitracy  parameters 

any  have  on  tJie  viscous  flow  near  the  stagnation  point  vlll  be  discussed 

la  the  latter  part  of  Chapter  H. 

As  esqpected)  the  esgpaasions  for  the  thexmodynsmle  properties 
exhibit  Crocco's  relation.  Cb  the  stagnation  streamline  the  flow  Is 
Iseatrppie;  Isentroplclty  is  nalatalhed  off  the  stagnation  streamline 
also,  vlth  the  exception  of  the  rotational  term  V  .  Ihe  pressure 
term  shows  no  effect  of  rotatlonsll'ty;  the  term  is  the 

centrifugal  pressure  gradient  due  to  curvature. 

To  complete  the  treatment  of  the  Invlscid  flcW)  the  velocities 
can  be  vrltten  down  by  substituting  (1  .§2)  into  (1 .4} 

u  =  Ax  [_i  +  ('^v-  + 


(1.23) 


CHumR  n 

Vliooua  glotr  ne>r  the  StMtnatlon  Point 
Oensnl  Coaalderatioas 

OSm  tesla  of  coosldazlag  th*  Tiseous  flew  iMur  ths  stacpmtion 
point  Is  Is^erstrom  and  Cole's  esqpanslOB  proesdure  (refsrenee  20),  as 
deserlted  briefly  In  the  Introduction.  Ow  essence  of  this  procedure 
Is  a  asgnlfleatlon  of  both  the  IndapsaAsat  coordiaate  and  the  velocity 
conponent  perpendicular  to  the  solid  surfhee  in  inverse  proportion  to 
a  "significant  viscous  length"  (l.e.,  the  bouBdaiyt>lBgrer  tMdoiess) 
and  a  velocity  based  on  It.  As  the  Unit  of  very  large  BqnMlds  nuniber 
.  Is  taken,  viscous  effects  will  be  Halted  to  a  vazy  thin  layer  near  the 
body  surface,  and  these  aagnlfleatlotts  then  pemlt  an  aaalytle  Investi¬ 
gation  of  the  structure  of  this  very  thin  layer.  Bils  layer  is  of 
courae  the  classical  boundary  layer  of  Prendtl;  and  lageratron  and 

(20) 

Cole's'  '  procedure  consists  of  laprovlag  the  bonadary-layer  result 
by  erpaadlng  all  flow  quantities  in  poeers  of  the  Inverse  sqnaze  root 
of  a  Bsynolds  number  based  on  sosw  slgalfloaat  length.  Bus, 
Ijteedlately  as  this  laprovsswnt  Is  eonaidarod  the  question  has  to  be 
raised  shat  this  significant  length  should  be.  Biera  is  no  such 
length  Inherent  in  the  elassioal  boiaidary>layer  solution  itself. 
Another  say  to  fozsulate  this  saaN  fne^tion  is:  shat  length  dhould 
the  boundary^layer  thickness  at  ths  stagaatlott  point  be  eonoarsA 
'  to  in  order  to  decide  shsther  or  not  a  eozveetion  texai  to  the  boundary^ 
layer  sblutlan  Is  necessary,  shea  eoMlAerlag  the  viseous  flair  near 
the  stagnation  point.  Ihia  is  a  ezusial  fusation,  shioh  sill  affect 
the  fonailation  of  the  entire  prOhlM* 
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OD0  Btloh  length  that  suggeats  Itself  la  the  radius  of 


curvature  at  the  nose.  Ibat  this  aagr  be  a  suitable  length  can  be 
seen  by  eonslderlag  a  stagnation  point  located  on  a  eurved  nose. 

Tor  very  large  Bsynolds  nuabers  the  boundsry^layer  thlekness  becomes 
negUglbly  small  cosqpared  to  the  nose  radius,  uhleh,  (by  geometric  Intu¬ 
ition)  Is  eq^lTalent  to  taMag  the  "Infinite"  radius  limit.  In  this 
Halt  than  the  claasleal  solution  of  flow  Inplnglng  on  an  "Infinite"  plate 
la  recovered.  On  the  other  hand  it  Is  easy  to  see  (relying  again  on 
intuition)  that,  as  the  Reynolds  nuibber  beccnes  mssUer,  and  hence  the 
boundaxy^layer  thlekness  larger  uhen  eonpared  to  the  nose  radius.  It 
any  become  necessary  to  consider  a  correction  to  the  claasleal  result 
due  to  the  curvature  of  the  boundary  lagrer  near  the  stagnation  point. 
(Subsefcent  analysis  will  shew  later  that  these  intuitive  considerations 
are  essentially  correct.)  If,  howrver,  the  nose  radius  Is  taken  a^  the 
sole  length  vhlch  could  be  of  significance  In  the  problem,  then  no 
correction  team  uhatsoever  can  be  adaltted  for  a  fla'U nosed  body* 

This  Is  certainly  contrary  to  expectation,  and  In  direct  contradiction 
to  the  presence  of  a  "boundary-layer  displacement  effect,"  as  dlseuseed 
In  the  Introduction.  Actually,  It  Is  not  necessary,  or  even  possible, 
to  decide  beforehand  shat  the  proper  reference  length  should  be.  It 
Is  clear  that  the  proposed  procedure  could  be  applied  with  the 
reference  length  left  axtltrazy,  and  ehatever  the  Important  reference 
length  of  lengths  may  be  they  will  appear  In  the  solution  upon  proper 
expansion  of  the  equations  of  motiom  and  proper  application  of  the 
boundary  conditions.  In  the  subsefoent  analysis,  the  nose  radius 
will  be  used  as  the  reference  length;  In  the  light  of  the  above  reamudu, 
this  eholoe  Is  one  of  coenrealonee  only  and  cannot  affeet  the  validity 
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of  tlw  outeone  of  the  aaelysls.  At  the  end  of  thie  elu^ter,  It  Hill 
he  poesihle  to  Identify  all  the  lon^aagmolds-mmiber  effects  In  texaw  of 
several  refezenee  lengths  to  vhleh  the  houndazy-layer  thickness  has 
to  he  coBvarad  vhen  considering  the  necessity  of  a  loi^Beynolds-nuaiber 
correction  to  the  houndary^lsyer  result. 


jDevelqpneBt  of  Theory 

Sie  procedure  for  investigating  the  viscous  flour  near  the 
stagnation  point  vill  then  eonalst  of  vrltlng  the  full  viscous 
eqiuations  of  aotlon>  as  derived  in  Appendix  A>  in  the  curvilinear 
houndazywieyer  coordinate  system  of  eqjoation  (l.l).  Nsgnifieation  of 
the  independent  coordinate  is  acccnpllahed  by  replacing  it 

with  the  houndary- layer  variable: 

y  =  '+1I  (2.» 

Magnification  of  velocity  VT  ,  and  the  desired 
ejqpanslon  of  the  velocities  is  aceonpllBhed  hy  vrltlng  the  strean 
function  in  terns  of  the  foUovlng  expansion: 


(2.2) 


Using  the  ^rasatry  of  the  prohlen  analogous  erpansions  can 
he  mitten  down  for  the  themodynaaic  properties: 


(2.3) 


p  =  ft  w)+'“] 

As  In  (l .6),  the  fluid  properties  are  assuaed  to  he 
teoperature  dependent,  and  can  be  esqjiaaded  in  a  Xsgrlof  series  about 
the  leading  tern  in  expansion  (2.3)  •  Ihe  dots  again  signify  deriva^- 
tives  with  respect  to  teaperature;  the  subscript  0  indicates 

evaluation  at  infinite  Reyixolds  number  and  x  »  0  ; 

I.  -  c..  t  ^.(T-T.K)  +  'i  CT-T,ll.f+,., 


2Q 

(2.4) 

(2.5) 


»  A,-*- 

»  ♦ 

t  =  V.,  +  i.,('r -'*i  ^  (T +  ••• 

(2.6) 

Using  ejqpeaslai  (2.3)  these  expressions  can  be  written: 


(2.7) 
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nie  mcnentuB  eolation  can  b«  vrittan  doim  "by 

eonalderiag  (1.1)  In  (A.7) 

>  f  5.  „  -1.  Hi: .  vts  + 

X  Vt  5jr  +  y.»y  40  ^ 

(«-sf  U»  i  (\*|7  R(^] 


('-trR"(— - 


-  /-  4L !«..  ^_ly  fc?i  -  - 


0**1)  u 
11(1*  l)  1 


+  -  tsl*-. ^  i*  +  - 


-v-*f ^ - 

'•'jjw  ^  ^  »«*>•  <>'^S  1A^>  1K>>-  ly 

■■r‘(.»*:i‘  R’t^  5  (T^ 

^  V  tTt  "■  V  ^  V  -^  '^'»^  ]  (2.8) 

Hoir  eiqpcuiBioiis  (2.2)  throuc^  (2.?)  and  (2.7)  can  1)e 
substituted  Into  this  equation,  and  the  results  grouped  according  to 
powers  of  snd  ^  Ihe  coefficients  of  the  R  ]fj  j 

constant,  R  »  *°d  teas  are 


,  and 


teas  aa 


S5.Ps  +K\:P.')=0 

SiPi  CV»?jp^  +  0 

ivpj  i-  v«v.'  ^  r.^f.  +  r»  V.  )'-i^  V*  *  0 
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b«nee: 


<f.  =  0 

V.  =0 


«■:  =0 


(2-9) 

(2.10) 

(2.11) 


—  V  =  - 

!.A'  *  V. 


(2.12) 

Qm  X  Bcnantum  equatim  li  obtalnad  by  using  (1*1) 


In  (A.8)j 


- ! _  /  S'  y  ^  ^  -  i-y  y  ^  - 

4.  K*  H/  -  1  J.  ^  %  - 

'  \  J  ^  ^  ^  ■“ 

_  _ I _  f  ft  s»  s  HUK  Sj^  _  ^  Sv 

0-5)^  5  3  0a-|)'-  ?  ' 

-  '*>  Slaiit  L*  .K  »  i  no«i>|.  Sjc  +  1  ^"-Ar  8v  , 

3  (Ttf^  S  iO»P*  ?  3  i  b  0-  iv  s 

+■4  -  -  '*'*'* ^  b  -  ~  _  y*>:>  Sx  4. 

M-l)  J  3  0^1^5  cTTi?  1 

+  yx>»»  _  Zr,  (^v  y^Ki.  ^  ^  _ 

v.^1  o  +  T'  ^  V  o  +  lr 

._  (uA>  j_  vK  I..  _  tu  ia  *>■ .  i  nf* /*-*  1.  y ««*»  _ 

—  ^  iCoT^  ^  ~U^ 

^  itiA  +  ^  u  + 

^  V-  -  ^T,,^  V  v'-=  -  <^>  ‘f  - 


(2.13) 


Bipeatlag  the  foxner  procedure,  eoeffldoata  of  the  X 


tern  give  the  foUoMlag  eqiietloa 


.  \ 


(2.14) 

Irene  ^  and  ^  in  the  above 

As  Ai 

eiq^reaaioci  are  fuaetloua  of  only,  the  exact  fuactlooal 

fom  dapeadli^  on  the  vlacoaity  va.  taaperatuze  lav  that  la  aaemed. 
a.  .«t  .Wtla,  1.  atlMd  Iff  cc^amt.  Of  tl»  2  J| 

ten  In  the  ^  aonentua  eqiuatlon  (2.13): 

+  ('n"-t)^,|i  +  Ov*<)M^  f!  -O+'^U.U  4!  + 


'  iJ  ^  «  \  \ 

2.  \»  -V  ^2.  =t 

^T,  S>n 


^  KU  ^  1^,  ]  ^  [1,“'  - 2 ^•+r-2  *2+.*  h  - 

-  b  *  2  is  ] + cn-i)  ^  I 


(2.15) 


2k 


Vtagjjy,  tlw  energy  equetloii  (A.11)  le  written  dom  In 
eoordlnatea  (1.1): 


R  C- 


tJi) 


iC'v.Pj-  'r,p.)]=0’iPi-f, 


>1) 


+  C.  ^  n""’  l,T,  +  0  *  D"  +  (H  ir  IT, 


+  O' 


ir'"  V.  j. 


I 


(2.16) 


The  dleelpetlon  function  ^  in  the  ahonre  expareealfla 
hae  not  been  expanded,  because,  to  the  order  of  terns  considered,  it 
vUl  not  contilbute  to  the  eguations.  In  order  to  ascertain  this, 
one  has  to  consider  the  form  of  the  dissipation  function  given  in 
(A.  10)  in  terns  of  the  general  coordinate  systen.  Sach  tern  vill  be 
of  the  fom  of  the  stream  function  occurring  twice,  and  four 
derlvstlTes  occurring  in  like  pairs.  In  order  to  contribute  to  the 
stagnatiooF>point  tern,  at  least  two  of  these  deriTmtlves  hanre  to  be 
)<  derivatives  (otherwise  the  tern  will  be  of  order  ^  ) . 

aince  is  of  order  *>nd  derivatives 

of  order  ,  the  largest  tern  in  the  dlasipatlcn  function 

will  be  of  order  viscosity  times  two  stream  functions  and  two 
derlvatlTss,  which  is  altogether  at  order  But  only 

terns  of  order  1  and  order  are  of  interest. 


How  (2.1)  through  (2.7)  can  be  substituted  into  energy 
equation  (2.l6),  and  terns  grouped  in  powers  and 

Ihe  leading  tern  is 


CoUaetlBg  the 


terns  gives 


i  'll 

H.  -  C'*«)(4.K  +  (*l  4.)  t^c...>U.4p;. 

Ti  ti'iSs 


(2.18) 


1. 


idiere  the  tezms  Y"  etc.,  are  s^caln  given,  bat  m  yet  on- 
specified  function  of  ^tl^)  .An  addltlcnsl  relatiaa  between 

the  variables  is  given  by  the  perfeot-gas  lav,  vhloh  was  assuasd 


■\?o  =  K 


-  44,'lri  + 

etc.  .X..  (2.19) 

low  the  fomulatlon  of  the  differential  eqiiatlons  is  ooagplete. 
Bqoatlons  (2.9)  throui^  (2.12),  (2.l4),  (2.1^),  sad  (2.17)  throu^^ 
(2.19)  define  two  coupled  sets  of  ordinary  differential  efoations.  In 
the  ser^th  and  first-order  variables  respeetively.  Om  serotbpoorder 
set  is  nonllnsari  it  is  essentially  the  boundary-layer  eguatioa  at  the 
stagnation  point.  33ie  first-order  set  is  Unaar,  with  the  eoeffioients 
and  ihhonoganeous  terns  ooaqposed  of  the  known  bouadazy>layer  solution. 
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In  ordar  to  proceed  vlth  tlie  solution  of  these  equations 
suitaUe  boundaiy  conditions  have  to  he  cpecifled.  Two  types  of  boundary 
coaditions  are  available;  at  the  asLid- fluid  interface  (i.e.,  the  "nail") 
and  in  the  Inviscid  flow,  vhere  the  viscous  solution  has  to  "merge"  into 
the  inviscid  solution  specified  earlier.  The  boundary  conditions  at 
the  wU  are  derived  In  Appendix  B,  and  are  repeated  here. 

'Veto')  -  0 
lot®'  =*  ® 


(2.20) 


and 


|,to)  ^  0 


>  a-  <  w  h-  I  ~  ^ 

'  ^  Hio) 

V  h"» 

(2.21) 

uhere  constants  K|  and  are  defined  by  the  above  expressions. 

In  order  to  detemine  the  boundary  conditions  in  the  inviscid  flow, 
one  can  use  "the  inner- flow"  independent  variable  (2.1)  in  the 
expressions  for  the  inviscid  flow  given  by  (1 .22)  throu^  (1 .2$) . 
Furthenaore,  as  explained  in  the  introducticn,  a  correction  will  be 
necessary  in  the  inviscid  ("outer")-flow  parameters,  due  to  the 
displacement  effect  of  the  boundary  layer  around  the  body.  OMs 
effect  will  change  the  velocity  gradient  A  ,  which  is  a 


27 


fuadaawntal  pamwter  of  the  pxoblom^'by’  an  usknom  Mwunt.  Thlt  ohuge 
can  he  ejqpxeaaed  aa  foUom: 

(2.22) 

Ohe  other  paraaeters  of  the  Isvlseld  flew  ulll  not  be 
ohanged  to  the  order  of  tezns  eonaidexed  here.  (A  aore  detailed 
dlaeuaalon  of  theae  eorrection  texma>  and  their  algnlfloaooe  win  be 
given  In  a  later  aection) .  The  invlaeld  atreaa  function  can  nov  be 
written  In  texma  of  the  vlaeoua  (i.e.,  "Inner”)  ewpanaion  aa 

T  =  iM  {7  *  iH  } 

^  (2.23) 

BlodLlarly,  the  expanslona  for  the  'UiemodjmaBd.e  propertlea 
In  the  Invlaold  flow  beccaw  ^ 

i  -  Ss.  I'  ^ 


P  =  ?s 


(2.2lf) 


Btpreaalona  (2.23)  ud  (2.24)  in  conjunction  with  (2.2) 
through  (2.3)  can  now  be  uaed  to  define  the  boundary  oondltlona  neeeaaary 
for  the  proper  "aergenee"  of  the  vlaeoua  "Inner"  flow  into  the  Invlaeld 
"outer"  flow.  Aa  beeoaea  very  large 


aB 


1 


1 


>  oo 


k(i)— ^  * 

2.9i. 


-1 


and 


— ►  0 

^  0 

-I?, (7))  — ♦  0 


(2.25) 


(2-26) 

"»  '7 

IMiig  tbeM  bouaduy  conditions  three  of  the  dlffexeBtlal 
equations,  nansly  (2.9)  throutd^  (2.11),  esn  be  integrated  iaaedlately 


k  ‘^1  =  0 


(2.27) 
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OwM  results  esa  Its  eoabiasd  with  the  pexleet>css  Isw 
(2.19)  to  oibtsln 


(2.3B) 


Qie  renalning  differential  equations  can  nov  be  sljq^fLed 
by  visiag  the  above  intemedlate  results.  Using  (2.27)  •ad  (2.2B) 
equation  (2.l4)  becones 


>» 


+ 


(2.29) 


The  X’  fflcnentuB  eqyiation  for  the  correction  texm,  (2.1§), 
is  slailarljr  oodlfled.  (trouplng  hoaogeneous  and  IxtbasogeneouB  texas  on 
the  left-  and.  rlgh't'haad  sides  respectlTely,  the  equation  beconws 


- 1 +.  ]  'H. 


^ 4»  4*4*  “ 

41'  .  1  Lti.  ^  V 

Ho 


^  (+: .  a  Vi  t .  V.  ^  'UK .  (tu:.  f\ ]  ^ 
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+  cif  ^1!  “  V*  J  ^  I  +  ■*“ 

+  v.  +Uo^*  •»• 

+  ^0  +a('-y.H',  ]  +  (vn-i)4*'  +M)U  ^  1 

H.  (2.30) 

nw  Btam  mibstitutioins  oan  Tm  aade  la  Uie  WMrgjr  «4imtiaas; 

(2.17)  iMCOMS 


\  C-«) \4\  ^  V  "  I  ^  ^ 

SlBllarly;  groiqplng  bonogtaisoas  aad  liduaogttMous  tezns 
(2.18)  'becooKis 


(2.31) 


(2.32) 

Sq^atlons  (2.29)  (2.31)  aw  eoutltut*  •  qr«t«a  of 

two  eooplod  BOB  linear  differential  equations  in  the  two  TazlaUes  4o 
and  ^  order  of  the  ooiddJMd  qrsten  is  S  »  thus  H 

boundary  eoodltlons  are  needed,  aquations  (2.20)  glTe  3  at  tbe  wall; 
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two  Bore  are  needed  “at  infinity.  “  Ote  proiper  IxMBdazy  eooditiooa 
are  then  froa  (2.25) 


f.tl) — *  I 

(2.33) 

Hieee  are  tlie  bcundaiy^leyer  eguationa;  their  solutlaae  can 
be  obtained  (nuaerleally)  as  socn  as  the  dependsase  of  fluid  properties 
upon  teaq^enture  is  specified.  OSie  wnlUtfr>«re*-stx«ni  stagtatlcop 
tesq^erature  ratio>  W  >  will  be  a  paxaaster  of  the  solution.  After 
these  results  are  obtained  the  solution  of  coupled  eqimtions  (2.30) 
and  (2.32)  can  be  considered.  Viaxlehle  %p,  in  equation  (2.30)  can 
be  eliMnated  by  Integrating  eqiiatlon  (2.12),  which  can  be  done  directly 
because,  using  (2.29) 


The  boundary  condition  for  >  as  given  in  (2.26),  can 

be  used  to  detenilne  the  miknown  constant  of  integration.  Int  the 
behavior  of  far  fron  the  wall  be  described  by 


T  ^  ~  -^.1,1  -» 7  -  7* 

where  is  detemined  fron  the  solution  of  the  boundarywiayer 

eqjuatlons.  Then,  far  fron  the  wall 


7”  i+*» 
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Brace,  ualag  (2.26) 


ead  the  expression 


(2.35) 


(2.36) 

estt  new  he  tised  in  (2.30). 

Before  proceeding  with  the  discussion  of  the  solution  of 
eqiustlons  (2.30)  and  (2.32),  result  (2.2B)  can  he  used  to  aodlfy  the 
houadsxy  conditions  at  the  nail  for  axid  as  given  In  (2.21) 

4,lo)  »  0 


(2.37) 

Ihe  tvo  constants  (ocr  rather  tmiflsl ration  of  constants)  lhat 

• 

cq^pear  In  the  shove  boundary  condltlcraaiw  azUtxazy  In  nngnltude;  their 
ratio  K,  /  Kj.  Is  also  azhltraxy,  depending  partly  on  SBQlvleally 
deteradned  soU^gas  Interaction  properties  ^  and  «<«w  .  Additional 

axtltraxy  constants  appear  la  the  houadaxy  conditions  for  "at 

Infinity"  as  given  In  (2.26),  anaely  aal  V  (this  latter 

only  for  the  \n  ^  I  case).  9w  linearity  of  tha  eqiratlons 
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■viggcvtf  tliat  It  Hill  iM  poitlbl*  to  ecnstruet  a  solutian  of  tba 
eqjimtloni  for  any  ecnbiaatloa  of  the  above  axUtrazy  eonstaBts  in 
tena  of  a  tm,  lAere  the  axbitiaxy  ecoataata  appear  aa  eoefflelenta. 
Baeh  of  the  fuoetloiia>  aaaoelated  with  aay  one  of  the  axbltraxy 
conataata,  oaa  be  solved  for  cooe  aad  for  all.  Ihe  aolutloBB  of  the 
eqpatloaa  eould  then  be  viltten  In  the  fora 

(2.3B) 

ahere  the  functions  with  subscript  c  are  associated  with  the 
solution  of  the  iniunogeneous  equation^  subject  to  the  boundary 
ccoOltlaoa  with  all  azbltzauy  constaata  vanishing;  whereas  the  other 
functions  represent  solutloos  of  the  hoanganeoua  eq|uatians>  subject 
to  boundary  conditions  associated  with  the  respective  eanstaats. 

Result  (2.36),  above,  then  Ixq^es  the  existence  of  ^  different  pairs 
of  functions  (4  pairs  for  n  =  0  ),  each  psdr  being  a  solution 

of  two  coupled  equations,  aad  subject  to  the  respective 

boundary  oondltions.  A  slapUflcatlon  Is  possible  If  one  d)serves  that 
the  pair  of  functions  aad  Identically  aatlafies  the 

hoBOgeaeous  eq^ations  (2.30)  aad  (2.32),  (because  the  left-hand  sides 
becQSK  Bcrely  the  derivatives  of  the  corresponding  b<Mndary^layer 
eqjoatlons  (2.29)  (2.312)*  functions,  nultlpUed 

by  aa  arbitrary  constant,  also  aatlsfles  the  slip  aad  tenperature-Jiaqp 
boundary  oondlUons  at  the  wall,  (2.|7)  prowidart  the  two  eonstanxe 
aad  are  equal.  !Ehis  solutloB  Is  coiviljaltfy  eqiulvaLan'  1 0 
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(23) 

XdJi  aai  'bouadazj^ljqrar  pcrtoitatton  telatlOK  tea  to  Blip, 

(24) 

or  tlia  altantata  aathod  uaad  ^  Muagltr'  '  of  "daprasBlag''  tlia  poaltloB 
of  the  BoU^fluld  Interface  In  the  houadary  lagrer  hy  a  length 
proportloBal  to  the  nean  free  path,  la  order  to  aeeount  for  the  Blip. 
looMller^^^  pointed  out  thB,t  depreBBlag  the  true  poBltloa  of  the 
nail  la  order  to  aeeount  for  the  allp.  Hill  not.  In  general,  aeeount 
eorreetly  for  the  effect  of  the  teeoerature  Juqp.  Caly  If  this 
"depreBBloa"  tiawe  the  local  teeq^ermtum  gradient  1b  eqnal  to  the 
teaperature  Jiap  alll  thla  q^proadoaticB  he  correct;  thla  oaae  1b 
equivalent  to  the  special  case  Ke.  In  the  teralnology 

of  the  present  analysle.  For  the  genaxal  ease,  > 

a  conreetion  function  can  then  he  detendaed,  nhleh  aseounts  for  the 
fact  that  the  two  eonstaats  are  not  eqnal.  IrpresBlons  (2.3B)  can 
then  he  rsgcltten  as  foUove 

^  ^  j;?  <•  ( VK.)  M!'"  i»]  V  c,  .  .M  iw 

(2.39) 

Bus,  to  detezadLne  the  affect  of  slip  and  teeQeratuze  Jap, 
only  one  pair  of  functions,  sad  ,  has  to  he 

deteaslaed,  instead  of  the  two  pslrs  la  (2.3B). 

Bw  differential  equations  and  houndazy  eondltlons  for  each 
of  the  functions  appearing  In  (2.39)  mb  nov  he  susHazisad.  First  of 
all,  the  goTsming  differential  efuatloBS  (2.30)  sad  (2.32)  earn  he 
aodlfied  and  nzltten  la  shortened  operator  notation  as  foUoers! 
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-I 


E  )  W  “  ^‘>•'*1^  (2.110) 

idiere  the  qnhbole  M  ead  E  denote  the  reqpeetlre  diffexentlal 
opevetore  for  the  naaentun  end  energy  egimtlooB;  end  and  E(. 
denote  the  reepectlve  Inhonogeneoufl  teae  iOTe>rl.ns  «  the  rid^t-hind 
■Idee  of  the  eQuatloQ.  Fron  (2.30)  aad  (2.32)>  the  above  operators  are: 


H  |M 
I  + 


*  i/'  (a..„ 


-h 


■nd 


e 

t  i-  Tt  +  V  ^  1  ^  ^  Z  r  ■* 

I,  V,  *^1  • 


+  +4f" 

V-* 


(2.1»2) 
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aim 


+  i  (\-Vy 


(2.43) 


+ [<''-''*A;l«-  *  s ^‘7*  <*•  -  ^  *!.  fe. 


_  3  b 


■ov  l«t  all  tte  aztltrmrjr  pKmwtan  In  (2.39)  vmnlah. 
9Mn  ttw  differential  eolations  to  lie  eolTed  are 


(2.M>) 


(2.45) 

Ibe  boundazy  coodltlcna  at  tlw  nail  are 

Vclol  -  0 

Vj«)  =  0 


\4-,cio)  =  0 


(2.46) 
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uhllCf  T«sy  ftiX  frcm.  Ham  mil 


1-*^  fa:  ^-o  4-1(11  -,- I 

for  Vl  1  4^4)  -»  0 

low  tba  corrmtloa  tezu  diw  to  tta*  ilip  and  taavemtare* 

JuQ  oftaetfl  e«a  be  eoneidered.  Bie  efuatiooe  are 

Ea.v,, 

(8.W) 

In  agreagaeat  mtb  (2.37)  Md  (2.39),  the  boondazy  eoodltlona 
at  the  mil  are 

=  I  (2.49) 

nie  "outer"  boundary  ecndltiona  are 

^ - >  oO 

V,  -*  0 

(2.50) 
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Diilag  (2.^3)  the  •qiiAtloas  for  tli*  tozai  aultlpllod  fegr  C>  are 


P  ^  J. 


(2.51) 

Sie  bouBduy  eoodltioao  at  tha  vail  axa  hOHOBtaeoua 


(^,.(01  =  0 

At  "Infinity",  using  (2.26) 


(2.52) 


^  oO 

h>  — >■ ' 


(2.53) 


Correction  functions  'V«>  and  can  ba  mitten 

down  in  tazBs  of  the  boundary>layar  aolutloiu  and  directly. 

To  obtain  these  eiqpressians  one  has  to  observe  that  eoaffielant 


appeared  in  the  invlscid  boundary  eondltlans  baeause  reloelty  gradient 
A  vas  ejqpandsd,  as  shorn  la  (2.22).  Oils  Texy  seas 

expansion  of  A  can  be  applied  to  tlM  ooigiaal  e^pxessioas  for 
the  strearihftinetlon  and  teaperatuxe,  as  given  In  (2.1)  throu^ 


(2.3)  •  Ibt  eonddarlag  fxaM  vltb  subsezlpt  I  ,  th«  cipaBaloa 
of  A  jritMs  ttaa  foUowlag  ucpMsolona 
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(2.54) 

OoMporlag  vith  (2.39),  cn*  ean  lonedlatoly  Idntiiy 

4»  =  iL4o  +  7V«*] 

=  i  1  (2-55) 

Substitution  Into  (2.51)  ttaxoui^  (2.53)  varlfles  those 
solutions.  Sxpresslons  (2.39)  esa  nov  be  written 

4.. ^  R [k,^.  ^  *'''4 

(2.56) 

FLaslly,  the  eqiuatlons  for  the  tens  proportional  to  V 


H(U,(‘w)=o 
E  4, V  4,4  =  0 


9m  bouadAxy  eoadltlons  at  tlM  mil  ax*  a^ala  hflaagaaaoua 

4iv  ^ 

|,v/  (o)  =  0 


(2.58) 


vhll*  far  from  the  vail 


Aj-#  -O 

ft^resalons  (2.U$)  through  (2.50)  and  (2.57)  throui^ 

(2.59)  give  the  differential  egjuatlaiM  and  the  aseoelated  houadazT’ 
eondltlonB  that  are  aeceseazy  to  deteznlae  the  five  pairs  of  functions 
(two  for  r\  m.  0  and  three  for  lA  I  )  that  oeeur  in 
eapressloiis  (2.56)  for  the  boundarTwiegrer  earreetion  texas  that  are 
under  eooslderatioii  in  the  present  ianrestigatloa.  9M  solution  of 
these  equations  is  contingent  upon  solution  of  the  stagnation- 
point  houndaxy-lagrer  equations  (2.29)  ■&!  (2. 31)  subject  to  houndaxy 
conditiQns  (2.20)  and  (2.33).  la  order  to  complete  formulation  of 
the  theoxy,  it  is  only  necessary  now  to  specify  the  dependence  of 
fluid  prqperties  upon  temperature.  This  dependence  can  now  be 
stipulated  as  follows  (cf .  Chapter  m) 

-<  T" 


i 


(2.60) 


becoM 


Btnee,  uslag  th«  ditflnltlons  at  (2.6) 


M-s. 


(2.61) 

Ihe  tflBipemture  dArlvatires  oeeurrlng  In  tha  aqiiuitlaoB  then 


L*  2.  Ii  ^  ^ 

f  =  00  ^0 


e^‘^.  (2.62) 

Ueiag  (2.6l)  end  (2.62)^  the  houndnijr-lBarer  equations  end 
the  opemtors  oeeurrlng  in  the  eqioatloni  for  the  eorreetion  terns 
esa  new  he  written  down  ejcpUeltlT'.  In  order  to  isolate  the  Id^iest- 


order  derivatives  the  expressions  will  all  he  divided  hgr  the  appropriate 
powers  of  .  Bouztdaxy-lagrer  equations  (2.29)  end  (2.31) 


heo 


HiL  4.- 

K'  c  V. 


III 


1  •* 


(2.63) 
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iriMM  unchMiged  boundary  eoadltlona  (2.20)  and  (2.33)  bar*  baan 


appandad  for  the  aaka  at  eoanraolent  rafareaea.  niaaa  aqpatloiis  ara 

eoaplataly  agulvalant  to  thoaa  originally  praaaatad  by  Brown 

and  than  aolvad  for  tha  twodlnenaional  eiaaa  by  Brown  and  Bonou^^ 

(4) 

Brown  and  Idvlngood'  ,  and  for  tha  aartally  ayonetrle  oaaa  by  Howa 
and  Ntraaan  (l6).  Tha  apparent  dlffarahoaa  ara  due  to  tha  alightly 
different  noxnallsatlona  of  tha  taaparatura  and  atraaai  funetlona. 

For  tha  correction  tazaw  <('•  ood  oparatora 

(2.1»1)  through  (2.44)  can  now  be  written  down 


e  ( i^, ,  -K.)  =  \  "  <*»+ 

+  ^(‘-0  7.]*.  ^  K'^  + 

+  di  It']  (a!.  +  (2-<^) 


EcU. ,  |i. v’-  ^  *  ^'- 


(2.65) 


niseusslQii  of  Sieoxy 

This  eaqpletes  the  presentation  of  all  the  differential 
equations  and  associated  boundary  conditions  that  axe  necessary  for 
obtaining  the  boundary- layer  and  first-order  boundary-layer 
correction  texais  for  stagnatloo-polnt  flow.  Ibe  ocnaeq]iisnces  of 
arbitrarily  choosing  nose  radius  R  as  the  referenee  length 
In  the  Seynolds  nunber  that  foms  the  basis  of  the  aqfagptotlc 


expansions  for  all  the  flow  and  thaxaodynanle  qiuantltles  In  the 

f«r 

problen,  can  new  be  discussed.  Let  the  viscous  lengthy  , 

utaleh  Is  Indicative  of  the  boundaxywiayer  thiekaess  at  the  stagnation 


points  be  denoted  by 


Oie  expansion  pameter  that  vas 

£ 

R 


used  in  the  above  analysis  Is  than  ;  the  ratio  of  the  boundary- 

layer  thickness  to  the  nose  radius.  As  R  becoMS  very  large 
conpared  to  ^  ,  the  first  of  the  correction  teas  In 

(2.36),  the  tea  ulthout  any  eoefflelent  and  labelled  by  subaexlpt  C 
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will  IMOOM  nagHgjlaly  aaall  eoagpartd  to  tbo  iMuiidaxT^ljgror  tom. 

Biis  tom  It  thon  ladood  duo  to  the  eurvmtuvo  of  ttao  noio;  it  irf^U 
not  bo  proooat  for  a  fla^noaod  body,  for  inotaaeo. 

I&  ordor  to  inveotigato  ttao  ptayoieal  olgnifleaneo  of  ttao 
renalalBg  eorroctloo  toms,  ttao  o^qioaoloa  paroaotor  ^  taas  to 
bo  eoosidorod  togother  witta  ttao  eoofflelonts  of  ttao  reapoetlTO  toxms. 

For  ttao  slip  and  toiaporature-JuBp  teiao,  ttao  paroaotor  ttaat  la 
slgBifleaat  is 


s  ar  irr  Ai  ^  ^ 

hi  s-  ^ 


vhore  ttao  constant  is  of  ordor  I  Sio  abovo  duors  that  this 

eorreetion  arises  vtaon  ttao  noon  froo  path  at  ttao  mil  boooaos 
signifleoat  ocoparod  to  ttao  boundazy^loyor  thieknoss.  Sineo  ^  boro 
is  in  ttao  denoalnator  ttao  botaarior  of  this  tom  is  guito  dlfforoat  fron 
the  fomor  one. 

Ibo  dlsploeeBoat  offoot  axisos  duo  to  a  etaaago  la  voloeity 
gradiont  "  at  saall  Beynolds  nuabors.  Itaoroforo,  ttao  length 

ttaat  Dill  bo  significant  in  this  Hayaolds  nuaibor  is  ttao  length  ttaat 

MAH 

plays  a  dotominiag  role  in  ostabllairiBg  the  aa^iltudo  of 
Oils  reforaaoo  length  L  could  bo  the  "also"  of  the  body 
pezpoadleular  to  tbs  free  straaa  for  the  ease  of  jttao  body  in  a  steady 
subsonic  straaa.  Or,  for  a  blunt  body  in  tayporaonle  flow,  it  is  actually 
the  nose  radius,  as  uill  bo  sbowB  ia  a  later  aoetion.  Oio  oorroction 
to  "A*  arises  ohea  ttao  boundaxyblayor  thieknoss  S  ia  largo 
onou^  eonparod  to  this  roforoaeo  length  L  TWmoan  eonstaat 

is  then  of  the  order  of  ttao  ratio  of  R  to  L 


t 


tbioh 


could  be  guite  Inalgnlfleant  if,  for  exanoOe,  tbe  stagnatioa  point 
is  locdted  on  a  "email  bump"  on  a  much  larger  body.  Qiia  shom  that 
the  curvature  and  displacement  effects  can  arise  Independently  from 
each  other. (See  sketch  below.) 


Esaaq^e  of  body 
with  dlspluemeot 
effect  negUglbe 
coosNured  to 
curvature  effect 
at  stagnation  point. 

finally,  the  last  tens  in  (2.^),  due  to  the  vortlelty 
effect  In  axially  syasMtrlc  flow,  can  be  codaaldexed.  It  Is  clear  from 
inspecting  the  coefficient  of  this  tern  that  xadlvm  R  cancels 
out;  Instead  ^  Is  compared  to  length  VV  .  A 
ihysleal  Intexpretatlon  of  this  length  can  be  obtained  from  the  linear 
velocity  profile  In  the  Inviscld  flow  given  In  (1 .23)«  as  Indicated 
In  die  sketch  below. 


Th*  ratio  of  tbs  bouadaxT^ligrtr  thlekMM  to  tbii  Iragth 
Mo) 

is  actually  Idantleal  to  Keap'a  '  vortlelty  paxMWtar;  the  ratio 
of  the  Tortieity  in  the  inviaeld  flow  to  the  "arerace"  Tortieity 
in  the  boundary  layer. 

These  coneideratiaaB  ahow  that  the  first-order  correction 
teas  to  the  stagnation-point  bowidaxy-layor  solution^  idiieh  are 
discussed  in  the  present  ianrestigation,  arise  due  to  four  different 
effects;  namely  the  curvature^  velocity  slip  and  teaoerature  Jwp« 
displacement,  and  vorticity  effects.  Each  of  these  effects  is 
associated  vith  a  different  loi^RBynolde-niaaMr  flow  parameter.  Three 
of  the  parameters  arise  from  a  comparison  of  the  boumdarywlayer 
thickness  with  three  different,  and  in  general,  independent, 
significant  lengths;  nose  radius  R  ,  length  U  (indleative  of 
body  size)  which  determines  velocity  gradient  ^  ,  and  length 

'/ V  which  is  associated  with  the  alppe  of  the  unlfom  inrlacid 

shear  flow  that  can  be  present  in  axially  syoBBetrie  stagnatiouppolnt 
flows.  The  fourth  paunmeter  is  eseentially  a  Xhudsen  number  based 
on  the  mean  free  path  neu  the  wall  and  the  boundary^layer  thickness. 

In  order  that  the  expansion  procedure  of  the  above  analyaia  be 
applicable  to  a  practical  problem  it  is  necessary  then  that  all 
four  Independent  parameters  described  above  be  an  order  of 
magnitude  smaller  than  one.  This  iaoUes  then  that  the  oadlest 
one  of  the  three  reference  lengths  be  comsidezably  larger  then  the 
boundary-layer  thickness.  It  also  iag^es  that  the  mean  free  path 
near  the  wall  be  considerably  maaller  than  the  boundarywiayer  thickness. 
This  last  reqnlreawnt  is  actually  implicit  in  the  use  of  the 


Hftvlaz^Stokea  equations,  because,  as  pointed  oat  in  the  introduction, 
these  eqjoatlonB  do  not  pexeiit  large  changes  in  flov  properties  (sueh 
as  occur  across  a  boundary  layer)  to  take  place  within  distances  of 
only  a  few  mesa  free  paths  long.  This  shows  that  the  use  of  the 
slip  axid  teoperature-jump  bouxidary  conditions  as  flrs'Uorder 
corrections  to  the  boundaxy^lsyer  eqirations  is  consistent  with  the 
use  of  the  Naviez^Stokes  equations.  If  ^  order  1  , 

the  ej^aasion  procedure  of  the  present  ansOysis  and  the  Ksvler- 
Stokes  equations  break  down  slnultaneously. 

If  all  four  of  the  ejqpension  parsaeters  are  sufficiently 
small  (i.e.,  smaller  than  of  order  I  ),  so  that  it  is  meaningful 
to  apply  the  Reynolds-xiuBibezwexpnaalon  prosednrs  deasxibed  in  the 
above  analysis,  it  is  interesting  to  eansider  wlait  type  of  terns  may 
arise  if  the  next  term  in  the  erpsnslon  is  considered.  The  energy 
and  manentum  equations  are  again  used,  as  in  (2.8)  throui^  (2.l8) 
now  terns  to  the 
collected.  On  the  left-hand  side  of  the  equations,  terns 
involving  die  new  (seconds order)  vmzidbles  >  Ve.  • 

will  appear.  Ohese  terns  will  all  be  linear,  with  the  coefficients 
composed  of  the  zero-th-order  (i.e.,  boundarywisyer)  terms.  Ihe 
linearity  of  the  equations  again  permits  splitting  the  seconid*order 
correction  terms  into  a  irjniber  of  separate,  mutually  independent 
effects.  These  effects  will  arise  due  to  various  parameters  that 
appear  in  the  Inhoangeneous  terms  and  in  the  boundsucy  conditions. 

The  inhomogeneous  terns  will  be  composed  of  various 
eonbinations  of  the  ssro^th  and  first-order  terns.  At  siost  two 
flrs'b'order  functions  can  be  ocntaLned  in  each  term.  8Lnee  the 


second  power  of  aqpansloa  variable 


H 


are 


1<8 


first-order  functlcos  can  have  either  of  three  axhltxmxT'  coefficients 
or  no  coefficient  at  all>  the  secondr  order  tezu  that  arise  ean  be 
associated  vlbh  conblnatlone  of  any  two  of  the  three  axbltraxy 
coefficients  or  no  coefficient  at  all/  the  secoodporder  terms  that 
arise  ean  be  associated  with  cfisiblnatlons  of  any  two  of  the  three 
arbitrary  coefficients/  or  any  one  of  the  coefficients/  or  no 
coefficient  at  all  (thereby  giving  rise  to  seven  different  effects 
for  the  axially  synnetrle  and  five  for  the  plane  flow  cases). 

Another  new  effect  will  also  appear  among  the  Inhomogeneous  terns  due 
to  the  presence  of  boundary-layer  terms  that  are  one  order  hlj^r  In 
the  expansion  about  the  x  =  o  point  than  the  stagnation- 

point  boundary-layer  terms  (l.e./  terms  like  ,  etc.).  ThuS/ 

to  this  second  order/  the  description  of  the  body  has  to  Include/  In 
addition  to  the  nose  radius  and  nose  wall  temperature/  also  the  rate 
of  change  of  these  quantities  with  In  order  to  specify  the 

stagnation-point  flow  (implying  the  appearance  of  two  more  arbitrary 
parameters).  As  mentioned  in  Cltapter  1/  additional  arbitrary 
parameters  win  undoubtedly  appesu:  In  the  "outer"  (l.e./  invlscidp 
flow)  boundary  conditions  to  the  second>order  correction  terms.  These 
conslderatlane  show  that  a  full  treatment  of  the  secondi-order 
correction  terms  would  be  very  laborftous  and  complicated  Indeed. 

!Ihese  correction  terms  win  not  be  considered  in  the  present 
analysis/  except  a  word  of  caution  has  to  be  added.  It  is  possible 
that  coefficients  of  soem  of  the  new  secondporder  terms  are  so  large 
that  the  aaqmaslco  procedure  of  the  present  analysis  (which  is  based 
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on  th»  noM  mdluo  u  the  slgoifleaat  length)  Is  not  appUoehle, 
even  thoun^  first-order  correction  terns  ere  all  sufficiently 
anali.  nils  iaplies  that  there  are  nsv  physical  paraoeters,  nhieh 
shov  up  only  in  the  second  (or  possibly  higher)-order  terns  in  an 
expansion  procedure  based  on  nose  radius,  but  vhich  are  of  the  seae 
order  as  the  boundary- layer  term  itself,  thereby  precluding  the 
possibility  of  this  lype  of  asyaptotic-ejqpansion  procedure.  For 
exaople,  the  rate  of  chsnge  of  curvature  at  the  stagnation  point 
could  be  so  large  that  this  effect  dcninates  the  viscous  flow  near 
the  stagnation  point,  even  thou|^  the  curvature  itself  is  not  large. 

For  those  eases  uben  one  of  the  expension  parsaeters  is  so 
large  that  the  eigpension  procedure  outlined  in  this  chapter  is  not 
applicable  any  more,  yet  the  Bavierw Stokes  equations  are  still 
applicable,  different  nethods  of  solution  have  to  be  considered. 

One  aethod,  vhich  is  feasible  under  certain  eireuastanees,  is 
solution  of  the  full  HavleivStokes  equations,  vlth  the  Bsynolds  nuniber 
based  on  radlxis  as  a  pareaeter.  Such  solutions  have  been  obtained 
for  hypersonic  flow  around  a  sphere,  with  a  constant-densl'ty  fluid 
and  concentric  shock  wave,  by  Frobstein  and  Kengp^^^^,  and  Hoshlzaki^^^^. 
Later  Hoshlzakl'  Included  the  slip  effect  for  this  case,  as  a 
separate  parsneter.  Another  aethod  is  solution  of  the  boundary^ 
layer  equations,  with  the  particular  strong  low^Bsynolds-nuaiber  effect 
as  a  paraaeter.  Rie  vortlci'^  effect  for  lacoopressible  fluid  is 

(19) 

considered  in  this  manner  by  Kenp'  .  A  soaswhat  different  approach, 
again  for  the  ^naaetrie  hypersonic- flow  case,  is  tuMd  by 

Oguchl^^^  and  Berrlng^^^^.  Their  aethod  consists  of  using  the  bow 
shock  wave  as  an  "outer"  boundary  condition  for  the  boundazy*layer 
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equations)  thereliy  oecountisg  for  lioth  the  vortlelty  and  dlaplaeaneat 

effects;  the  foxner  for  a  S  ^  eoosteat,  the  latter  for  a  /.j  =- 

(29) 

ooostaut  fluid.  More  recently  Oguehl^  *  presented  an  aaalytie  solu¬ 
tion  for  this  ssoe  axially  synaetrlc  constant- density  hypersonic- flow 
case>  in  tents  of  an  expansion  in  shock  density  ratio  e 
This  analysis  includes  the  curvature,  displacesient,  and  vortieity 
effects  siimiltaaeously.  Finally,  the  method  of  Frobstein  and 

Ksnp^^^ )  also  been  applied  to  a  variable-property  fliild  by 
(32) 

Prcibsteln  and  Ho  again,  for  the  case  of  a  sphere  in  hypersmic 
flow. 

In  the  present  analysis  the  restriction  to  small  values  of 

the  panoeters  vas  accepted,  as  this  disadvantage  was  eoapensated  by 

the  possibility  of  identifying  the  four  first-order  lowkRsynolds- 

nundber  effects  ^  cooparing  them  on  an  equal  footing.  The  validity 

of  the  present  approach  can  be  extended  by  construetlng  "hybrid"  or 

"coaqtosite"  solutions,  including  only  the  largest  of  the  above 

effects  in  the  nonllpeeLr  solution,  and  accounting  for  the  remaining 

loia>Iteynold8-number  effects  by  the  pertuxbatlon  procedure  of  the 

(l9) 

present  analysis.  For  exanple,  Keop's  '  parametric  solution  of  the 
boundaiy-layar  equations  subject  to  a  vortical  outer  boundsuy 
condition  could  be  "perturbed"  vlth  respect  to  slip  and  or  curvature; 
etc.  Such  a  solution  could  be  useful  in  a  flow  regime  vith  large 
vortieity  interaction  effect  but  coqpau»tively  small  slip  and 
curvature  panoaaters.  The  most  useful  approach  mill  be  different  for 
each  problem,  and  can  be  determined  by  estimating  the  orders  of 
magnitude  of  the  various  low^Rsynolds-nuMber  flow  paraasters. 


In  order  to  apply  the  results  of  the  foregoing  theory  to 
flows  around  specific  bodies,  the  flow  parameters  that  are  Inherent 
In  the  stagnation-point  flow  problem  have  to  be  determined  for  the 
specific  flow  exaaqple.  Some  simple  examples  and  applications  will 
be  considered  In  the  following  paragraphs.  In  some  cases,  it  may  be 
more  appropriate  to  use  experimental  results  to  determine  the 
parameters  under  consideration. 

Subsonic  flow  around  circular  cylinder  and  sphere. 

At  low  Mach  numbers,  the  results  of  Incompressible  flow 
curound  these  bodies  can  be  used.  Ohe  stream  functions  are  given  for 
the  cylinder 

^  (3.1) 

and  for  the  sphere 

Using  the  boundary* layer  coordinate  system  of  (l  .1)  and 
expanding  about  the  stagnation  point,  these  expressions  becosM 
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(3.2) 
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Oovftxliig  with  (1>22), 

(3.3*) 


V  =  0 

(3.3b) 

Do*  to  the  presenee  of  tiuctoulenea  and  tha  poaaibllltgr  of 

aapumtloiif  dlqplaeaMat  eonstaat  hu  to  ba  laft 

luidatazsiaad  for  this  tgrpa  of  flow,  lha  r—alnlng  pamatara  axa 

asaantlally  axtltxmzy  thamodyimle  propartlaa. 

^yparsonle  flew  around  eireular  eyUadar  and  ipbara. 

Invlaeid  solutlona  for  bgrparaonle  flow  arouad  eireular 

(4-3) 

eyllxider  and  aptaera  uara  givan  bgr  TULtlaua^  '  and  alao  Bagraa  aad 
Froibataln^^^  for  the  cyllxidar,  and  Idi^tblU^^^  for  the  apbara.  Both 
of  thaaa  aolutlooa  are  predicated  upon  three  aaauavtlona;  ahoak  abi^ 
that  la  eonoantrle  with  tha  bodgr»  laeoaoraanlbla  fluid  behind  tha 
ahoek,  danal'ty  ratio  aeroaa  shock  la  a  eonstaat.  Xha  stream  functions 
are,  for  the  ejrllndar 

(3.11.) 

which  daflaas  and  ;  and  for  tha  «phez« 


for  tha  eylindar  aad  sphere  xagpaetlTaly. 
for  tha  eyllnder 


aad  for  the  sphere 
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\t  »««.  V  I  \ot.  K^.  J 


(3.4b) 


vbere  «  It  the  density  ratio  aerost  the  daoek  ^  > 

and  the  radius  of  the  hhoek.  Let  the  Imrlseld  Shook 

standoff  dlstaaee  be  4  ,  so  that 


(3.5) 

In  the  boundary^ layer  coordinate  systan  of  (1.1)  then, 
and  eonslderlag  the  leading  power  In  x*  only,  the  expressions 
for  the  respeotlTe  straaa  funetlons  beeona,  for  the  cylinder: 


=  e  U, 


f  = 


and  for  the  sidiere: 


(3.6a) 


^  0-e)(>-te) 
\5’e. 


(3.6h) 

Since  the  body  outline  Is  part  of  the  stagnation  streanline, 
the  expressions  In  brackets  In  (3.6)  are  equal  to  0  at  0 

a 

Ibis  relates  to  e 


(3.7a) 


3  _  5  C\-  <vt)  ^  iCv-OC'- i) -■  0 


(3.7b) 


for 


IhAM  laq^Uett  r«latlcn«  ein  bt  «qpaad*A  in  a  Mxlat  toot 
e  (ef .  Buyas  and  Frobttaln^^^ 


(3.8a) 


(t  r  tf.-.  ] 

(3.8l>) 

How  stream  functions  (3*6)  can  be  expanded  in  powers  of  ; 


't  «[■»  - 

Agaln«  eosqparlag  wltb  (1.22),  tbe  velocity  gradient  for 
tbe  ejrllnder  becomes 


(3.9a) 


(3.9b) 


Similarly  for  tbe  ^bere 


(3.10a) 


V  ^  i 


(3.10b) 
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In  this  case,  the  displacement  effect  can  also  he  ascertained.; 
the  shock  standroff  distance  will  he  changed  from  the  Invlscld  value  due 
to  the  presence  of  the  boundary  layer^  without  (to  this  order) 
affecting  the  shape  of  the  shock  surface.  It  will  then  appear  as  If 
the  nose  radius  were  Increased  hy  the  houndary-layer  displacement 
thickness  ^  .  Then,  In  the  velocity  gradients  In  (3<10) 

above 


1 


Comparing  with  (2.22),  it  Is  apparent  that 


(3.11) 


the  corresponding  low>Beyn61dB-nimd)er  correction  to  the  shock 
stand-off  distance  is  also  of  Interest 


For  small  e 
leading  tenis  are 


(3.12) 


esqpiesBlons  (3.10)  can  he  ejqpcuided;  the 


(3.13a) 

for  the  cylinder,  and 

A  .  +  CTi.’'*!  .] 

V-  Ufi** 


(3.l3h) 
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for  th«  sptaoro  z«ipeetiv«ly.  dimplif  of  ■tud'Off  ilartsMwt  , 

voloeltar  gradloat  ,  and  vortiel'^  panuMtor  V  art 

ahoMB  la  Ilgurea  1  aad  &,  for  ttaa  aaka  of  eonranlant  rafaranea.  Sw 
oorraapoodlag  xuoiariaal  eoiiputatlonB  ara  ahom  In  AppanftUc  C. 

Bapirleally  datemlned  stagnation  point  valoeHy  gradlant 

va.  fra^atraan  Maohpsuoiber  eurvas  art  glvan  "by  Badiotko  aad  Backaltlx'  ^ 

l3k) 

for  elreular  cyllndars,  and  Oy  Crawford  aad  MeCauley'  ' ,  aad  BcacLg'  ' 
for  spheraa. 

nig^t  of  a  SLunt  Body  throagih  tha  Aiaoaphere. 

An  espaelalXy  Intaraatlng  aad  praetleaUy  algnlfleaat 
^KpUcatlon  of  loi^B^mOldB-nuBibar  atagnatlcu^polat  flow  oeeura  In  the 
eaae  of  a  iMdy  flying  at  hl^  altltudas.  In  order  to  datandna  wbathar 
any  of  the  four  correction  teas  daaerlbed  In  tha  analysis  of  Chapter  H 
are  nacassazy  or  i^rpprlata,  one  has  to  aatlaate  tha  magnitude  of  the 
four  req^tlva  paxamstara  at  various  fllid^t  apaads  aad  altitudes  for 
a  body  of  given  size.  In  order  to  detamljie  the  houadauywiayer 
thlekness  at  the  stagnation  point  tha  stagnation  kinaantle  vlaeosl'^ 
aad  velool'^  gradlant  A  have  to  ha  knomn.  Om  fozawr  can  be 
datezalned  unequivocally  for  a  given  altitude  aad  flii^t  spaed. 

Velocity  gradlant  A  will  ha  proportional  to  »  where  L 

la  soma  significant  lengthy  usually  ladleativa  of  body  size.  Leaving 

this  body  length  uzwpedfled,  boundary«layar  thlclmaaa  will  be 

17, 

proportional  to  qjuaatlty  ,  which  la  a  funetion  of  altitude 


aad  fU^it  spaed  only.  In  oxdar  to  detazmiiie  wbathar  tha  eurvatura 
aad  dlq^laoeaant  corraotlon  affects  ara  laportaat,  the  quantity 
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tbatlji  b*  divided  bgr  the  •q^exe  root  at  the  leiveetlve  lei^th* 

the  two  Icagthe  "r"  ead  "L"  tarn  at  iifferest  oedere  at  aignltude 

Vu  rf, 

ie  the  proper  peit— ter  to  detexMlm  the  etgnifleeBee 
of  the  eurmtaire  effect,  vhloh  is  then  dlffemt  free  the  dlcpleeeBeBt- 
effeot  pexeaeter 

To  aeeees  the  elgBlfleanee  of  the  ilip  Md  tenpemture-Jtnq? 
effeetc,  the  neea  free  path  at  the  alaiffvtlGK  pedait  IM  to  be  eenpared 


to  the  bouadary^legrer  thiekneee.  The  foaewr  een  be  «Qaulated  from 
the  thezieodarBaBle  propeorties,  ead  thus  sen  be  detefldJMd  as  a  fuaetlon 
of  ^  and  H  ooly.  Tfaue  the  ratio 

L‘  -  -ii- 

i  Jf  (L 

d(q^eiids  OB  a  paxaeeter,  uhleh  Is  a  fuaetlon  of  V  aad  ^  only, 
divided  by  the  square  3root  of  referenee  length  \_  (TOr  a  strongly 
eooled  body,  the  anaa  free  path  at  the  wU  omy  be  eoneldezably  saadler 
than  the  aean  free  path  at  the  stagflation  condition,  the  effect  of  this 
will  be  discussed  la  detail  la  Gtaeptar  17). 

For  axially  synfletxle  blunt  bodies  la  supersonle  flow, 
the  vortlolty  correction  effect  aay  here  to  be  considered  due  to  the 
preaeaee  of  a. curved  shock  mve.  At  hl|d^  Nfleh  audbers,  Idi^ithlll's 
^pproodnatlon  (discussed  la  the  previous  seetlon)  Is  reasontsbly  accurate 
If  the  aoae  outline  aad  the  ahoflk  are  oonceatzle.  If  this  la  aot  the 
oaae,  an  enplrleally  (or  othemiae)  detexnlned  ahoek  radius  oaa  be 
used  Intttad  of  the  nose  radius,  aquation  (3*10)  flhoMi,  that  la 
addition  to  zadlua  R  ,  vortlelly  pamaater  V  la  also 
dependent  on  shock  density  ratio  e  ,  iMah  la  a  fuaetlon 
of  ^  aad  H  only.  Thus  the  paraaeter  that  deteadaea  the 
aa^dtttte  of  the  vortlolty  effect 


!■  thn  tcMplef  ly  depodant  on  V  and  H  vlth  tte  «»o«ptlon  of 
a  fOetor  of 

Omso  eonaldaratloiui  dior  that  It  la  poaaiUa  to  plot  on 
a  “V  T8.  H  ehart  fanUlaa  of  liaaa  ahowlog  tte  aacaltuda 
of  all  four  of  tte  aigDifleant  pamatars,  laanrlag  out  tte  affaet  of 
body  alM  aa  a  aultlplleatlTa  faator  of  aztltrazy  nagiltuita.  Sudi  a 
plot  la  ataott  la  llgura  3«  Ite  foUoalag  finlllaa  of  Uaaa  ara 
hhoaa 

i  j  “rt 

&  t  tterafora  \  aad  V 

r?,  (fi 

vrU  ,  P 

i  I 


(3.)k) 

Bia  plot  la  baaad  on  tte  AHDC  nodal  atnoaptero  (rafaranea  42) 
raal-gaa  affaeta  ara  laeludad,  la  tte  ealeulatlona.  Ite  proeadora  aad 
maMxtaal  datalla  ara  glraa  la  Appandlx  S. 

Jar  a  glraa  alaa  body,  tte  Uaaa  of  figura  3  aaa  ba  uaad 
to  dallaaata  tte  raglon  of  ivpllaaUllty  of  tte  axpaaalon  proeadura 
of  tte  prasaat  aaaJyala.  All  four  of  tte  axpaaalon  pamaatara  hara 
to  ba  laaa  than  1 .0  (aetuaUy  unlaaa  ttey  ara  laaa  ttea  about  0.3« 
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th*  trror  eoaiitted  in  neglaotlag  tbn  ■Moii(WQrter  t«BM  Mgr  b*  aors 
tbaa  10)1).  laapeetion  of  th*  flcuxa  Mom  tiatb  &t  lav  RAvonle 
speodv  tlM  eurvnturo  (nd  dlypOLn— t)  eoorvMtioa  •ffveti 
beeoM  sigaifleaat  at  a  ouch  lowr  altltudv  tban  liw  slip  «ad 
taqwzatuza- jtOBp  effects.  On  the  otlMr  taaiid>  at  hgrparsonie  vpMds 
the  slip  eorreotlon  occurs  at  a  s&ssidiat  lower  altitude  than  the 
evirrature  sad  dlsplaceoent  correetions.  Tat  vptaeres  at  hypersonic 
speeds  the  vortlcity  correction  effect  is  largter  tor  orders  of 
nsgnitude  than  the  other  tarns,  eqpeeinllgr  as  the  abode  dsnaitgr 
ratio  hecomee  onall.  Also,  id  appaara,  that  as  the  hnuartaTy»layer 
sad  ahocl^lagrer  thicknesses  heeone  of  the  sane  ozder  of  aagnitude, 
the  altitude  is  already  too  hidh  (Md  hsooa  the  Boraolds  nitber  is 
too  low)  for  the  appllcahiUty  of  tfta  eaqpMaioB  prooioAse.  A  very 
detailed  and  thorou£^  dlsousaioa  of  Ite  sueesaalve  flow  ragiaas 
at  hypersonic  flow  has  been  given  recently  ky  Frobstain  and  ^ ; 

there  is  no  need  to  reiterate  these  results  here  exoept  to  point  out 
that  the  paraaeters  plotted  on  figure  3  eleaxGy  iadleate  this 
succession.  The  expansion  proeeduze  of  the  present  annlgrsis  is 
applicable  near  the  higtwBeynnl di^mafter  (i.e.,  low>altitude)  end 
of  this  spectrun;  especially  for  the  axially  ^naaBtrie  ease,  etwre 
the  vortlcity  "correction"  hecones  as  large  as  the  boundazy*lagrer  tern 
Itself  at  conparatively  low  altitudes,  Mere  the  other  corrections  are 
still  aaall.  At  hypersonic  flow  then  either  the  "eaaat"  eolations 
for  a  sphere  by  Freibstein  and  or  by  EoMLsaki^^^^  dhould  be 

used  (because  the  vortlcity  peuceneter  is  too  large  for  the  applicability 
of  the  present  eigpanslon  procedure);  or,  at  low  enoo^  altitudes. 
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idMx*  tbe  ttpusion  proeadur*  la  ■ppOloablti  only  th*  rortleitgr 
eosTBetlca  Is  slgolfleant,  sad  therefars  tbs  'hHdLflsd  bouBdaxy>lagrsr 
solvtlaBs''  (idMxa  the  rortieitgr  bouadaxy  eondltloB  is  a  paxsaster)  of 
C|guehl^^^>  and  Bezrlxig^^^^  are  equally  praetieal  (this  Is 
the  "Tortlolt^lntexactloa  regjLne"  of  Frobateln  and  Ilffire  3 

ladleates  that  the  flov  x«gloa  vheze  all  four  of  the  eorreetion 
effects  are  about  equally  large  oeeura  In  the  supersooie  regins 
(betwihahout  2000  and  6000  ft/see).  Zt  la  on  this  basis  that  the 
fluid  properties  used  In  the  solntiA  cf  the  equations^  as  presented 
In  the  next  chapter^  vere  chosen.  Die  proeedure  for  obtaining  these 
properties  Is  Shoim  In  Appendix  S.  Oia  bo^r  else  irfU  iddft  Ihs  flew 
regtsas  on  the  altitude  seals,  but  not  on  the  velodly  scale,  sinee 
all  panaaeters  are  divided  by  the  s^axe  root  of  the  al0alfieaat  length. 


caitfgn  If 


m—gj— 1  BBsuitfl 

n«sentatloix  of  Results 

OSis  dlffsrentisl  eqjufttioBS  prsaeated  in  Cteptsr  H,  idth 
fluid  properties  given  In  Appendix  K,  vex*  solred  masexlaallgr  at  tbe 
Comill  (Xaq^tlng  Center  by  wssns  of  a  Burrora  220  eleetronle  oo^puter. 
Soae  details  of  the  cooputlng  proeedure  are  given  in  Appendix  0.  Bia 
results  of  the  eonputatlons  are  deteznlaatlonB  of  "serott^"  and 
firs-Uorder  terns  (l<e.y  boundazy^layar  aad.  flrat^ooedar  *aane|iEl.dn. terns) 
for  the  noordlnetisloaal  strain  function  ^  and  tn^pexatnra  function 
>j|4  .  In  addition  to  the  functions,  ttadr  derivatives,  up  to  an 
Including  the  highest  that  occur  In  the  differential  equations  (i>e., 
the  first  three  for  4  *ad  the  first  two  for  ^),  ware  also  eonputed. 
Tiiw  obtained  solutions  vere  •*•->»»"  usad  to  dstemlna  sons  additional 
quantities  of  practical  interest.  Itasse  quantities  are  the  two 
veloeltgr  eoeponents,  the  two  aasB-velodty  eonpcnants,  tegq^rature 
end  density  profiles,  and  the  variation  of  vortlcity,  diear  (parallel 
to  the  wall),  aad  heat-tnnsfsr  rata  (aoranl  to  the  wall)  in  the 
TlseouB  layer. 

RaealUng  the  deflnitlonB  given  In  (2.2)  and  (2.3),  one 

can  write 


JV _  =  V 


1  _  T.  ,  i  I 

■'i 
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vUeli  AiflMS  >  Vo  >  %  >  "^0  * 

M^TCiiion  eaa  b«  nrlttaa  for  tl»  dnsit7>  ‘hy  usiag  rtwlt 

(2.38)  in  {2.k) 


<i  S.  ^  i  1> 
S,  '  i  S, 


'  -  i 


(1^.2) 


ATI  th*  otbAr  qjiaiititl**  of  intorost  on  Tio  wlttoa  domi 
In  •»*>»  MOMr  hr  "■<"8  ^  «qp*oi«loM  lad  ta»  rolntiaui 

doTClOPod  in  iBpMdleoo  A  and  B  BotuMn  tlM  •trmm  ftneUoa  aiiA  lii* 
Tuloiui  flow  foantitle*.  Baa 

V  .l-i-  .1\4S 


’1^-1 

.  s"  ^  «5i?  .Iff  C5!i'  ^  V» 

-TT/^  ^ 

^  ^ \*i  ^  ^  *  -  ’^’1 U* 

iM  -  ^  =  -  -v  I  =  V*  (4.3) 
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A  furtbar  braakdoim  of  Ite  iImt*  MgacotsloBB  ooonrs  idMii 
tbe  dlffortnt  effaeta  ahlch  aaka  iqp  tho  llrat»OKdar  eoaraetlai  tan 
(as  darivad  in  Chapter  n)  ava  eonalAanA  aapaxmtaljr*  To  faailltate 
this  prooeaa,  let  the  tezsia  due  to  veloeltgr  allp  and  taqparature  ^uap 
he  radsflaed  as  foILoira 


(4.^) 

(It  ngr  he  noted  that  the  ten  aaxtod  Iqr  shbeerlpt  >L  la  due  to 
both  the  veloeltgr- allp  and  taopenture-Jav  houBdaxy  eondltloaa,  and 
the  tan  oaited  by  J  arlaea  not  due  to  the  taapentora-Jiap 
houadaxy  condition  alone,  hut  nther  due  to  the  dlfTannee  hetwsen  the 
oonstanta  of  proportionality  for  the  valoel'^allp  and  tiewpenture- 
Jnp  houndaiy  oondltlons  to  the  raapaativa  sndlanta  at  the  vail) . 
Ualjv  this  definition  and  the  raaults  of  Chapter  n  the  straan 
function  ean  he  expanded  aa 


(^^.5) 

All  the  quantltlas  In  (4.3)  ean  he  vrltten  doan  In  like  aaimer.  S.g. 


0  .  'i*  +  - 


.  1  + 


Vo  V»4L  ^  C 1?; '  ‘^  ^  '^'* 


4-  v'V  ^^oVw  +  \o  4AnV 


(1^.6) 

etc. 

Die  results  of  the  nuaexleal  ealeulatlons  are  given  In 
Tables  I  eal  U.  liable  I  gives  this  ”lneo«presslble''  (i.e.,  oonetent- 
Qai^property)  results;  the  boundexy^iayar  texa  (given  for  reference 
only)«  curvature- eorreetion  ten,  dlsplaeeawnt-eoxraetian  ten, 
velocity^  slip  end  teapezatur^Jnp  ten,  and  the  voartldtyweorxeetion 
ten  are  tabulated  In  that  order.  Tor  each  ten,  the  two-dlaensicnal 
(vA  =.  0)  and  axially  qnawtzle  (n>  i  )  results  are  tabulated 
sub8eq^ently,  except  for  the  vortldty  correction  ten  vblch  exists 
only  for  the  Mri«^ny  sjnaetzlc  ease.  Oie  appropriate  differential 
equations,  boundary  conditions,  and  fonulae  e^presalag  all  the 
quantities  listed  In  (^.3)  axe  given  as  a  eonvedent  reference  at 
the  appropriate  sections  of  the  table. 

A  tabulation  for  the  "eonpreadble"  (l.e.,  varlable- 

flul^property)  case  Is  given  In  Ibble  n.  Basults  are  shown  for 
tiM  vV  ■=  0.73,  0.3,  0.23,  and  0.1  oases,  flia  W  =  1.0  ease 
Is  Identical  to  the  "Inoonpreadble"  result  (ef.  Appendix  7).  To 
iUustnte  these  results,  velodty  end  tenpenture  profiles  hare  been 
eoaputad  and  axe  plotted  In  llgaxes  4  and  3*  £Bioan  are  the  uneorreeted 
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teunduT^ljqrvr  profiles^  bouadaxT^liartr  profUts  ciorraetad  dua  to  tha 
effaet  of  eurratuza  only^  diaplaciwit  affaet  0Bly>  raloeity  allp  and 
toapezatojra  only,  aad  rortlei'^  affaet  ooly.  Far  aaeh  eaaa,  tha 
partiita&t  axpaaalon  paruNtar  aas  aatnaad  to  ba  a.g.,  for  the  euratp 
tura  eorraetlon  only  profile  [J  z  Uf,-h  0.2  Uic 

For  the  yaloelty  slip  aad  taqperatura  Jiap  eorraetlon  taxa  the  pareaetar 
aas  asatawd  to  be  1.0,  idileh  is  tba  rti^t  order  of  aaipiltuda.  To  aaphaslse 
tbe  affaet  of  ecaprassiblll'^,  e^bllag  ratio  W~  OJ  «as  dtosan  la 
these  aaaplas.  Tbe  obtained  eurvas  ara  shoan  la  tao  sapazmta  grouplags. 

In  figures  ta  aad  5a,  the  veloel'ty  aad  temperature  profiles  respeetlTaiy 
ara  grouped  together  according  to  the  type  of  eorraetlon  that  Is  (or  Is 
not)  considered.  Buis,  separata  groups  are  glTea  for  the  boiiadary>laarar 
profiles  only,  for  bouadary^layar  profiles  eorraetad  for  curvature  only, 
for  dlsplaceamit  effect  only,  etc.  Baeh  group  of  four  eweres  then  brings 
out  the  difference  betaeen  the  two  dlneaslonal  aad  axially  syaaetrte  pro* 
files  for  both  the  lacoavresslble  aad  eoapresaible  (i.e.,  constant*  aad 
Tazlable>fluldr*property’)  cases.  In  the  second  grouplag  (Figures  4b  aad 
5b),  the  profiles  are  grouped  together  according  to  the  type  of  flow,  i.e., 
all  two*dlaenslooal  iaccapressible  profiles  are  ahowa  la  the  aeae  group, 
etc.  This  grouping  then  brings  out  tbe  differences  between  the  various 
eorreetlon  effects,  conpares  then  to  each  other  and  the  uneorrected  (i.e., 
bouadaxywXKgrer)  profiles.  The  large  increases  la  fluid  velocity  aad 
taaperature  at  the  wall  due  to  the  veloelty^sllp  aad  teaperature-junp 
boundary  conditions  are  especially  noteworthy  for  the  ccnprassible  case. 
This  large  increase  la  due  anlaly  to  the  terns  with  the  subscript  la 
egoation  (4.4),  i.e.,  to  the  terns  arising  ftron  the  difference  between 
constaats  of  proportionality  aad  .  The  resulting  velooity 

aad  tngperature  valnsa  at  the  wall  are  (aeeordlag  to  the  eurves  of 
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figures  k,  S)  naarly  eq)UAl  to,  or  BUgbtly  excooftlag  tbe  freo-atroui 
values;  also,  the  profile  slopes  at  tbe  nail  appear  to  reverse  sign. 
Such  a  large  ebaage  cannot  be  esqpaoted  to  oeeur  la  tte  actual  flow; 

It  shows  that  for  the  parameters  that  were  assuesd  v<  =  o-i  >  'i'-i  >  \ 
the  expansion  procedure  of  the  present  analysis  cannot  he  used  any 
more  If  the  velocity-slip  and  tea^ratur^Juap  expansion  parameter 
\  N ,  Is  20it  or  larger. 

qm  qiuaatlties  that  are  of  tbe  most  praetlaal  interest  are 
the  shear  and  heat>  transfer  rate  at  the  sail  (  ^  )  Table 

III  gives  these  quantities  for  the  range  at  eooUag  ratios  that  «er« 
eagployed  in  the  caleulations.  Tbe  hea'Utraasfer  rates  are  expected 
to  be  very  much  dependent  on  the  temperature  dlfltaanee  between  the 
free  stream  and  the  wall ,  *  s '  ^  *  A  direct 

prcqiiartlooality  of  wall  heat>traasfer  rate  to  this  teaperature 
difference  is  usually  assumed;  then,  in  order  to  properly  naasallze 
tbe  noopdlnensioeal  wall  heat-transfer  rates,  they  hhould  be 
divided  by  the  quantity  \  -  W  The  haat-traaafer  rates  are 

therefore  given  in  this  aomnllsed  fooa;  they  are  egcdvalent  to  tbe 
usual  b<iundazyblaye]>heat-traasfer  parameter  ^  where 

is  the  Russelt  nuniber.  A  plot  of  shear  and  heat-transfer 
rate  versus  cooling  ratio  W  is  given  in  Flgares  6  and  7 
respectively.  For  the  boundary-layar  term,  it  is  wall  knosn  (e.g. 
Lees'  ')  that  the  sheaur  is  not  very  sensitive  to  the  presence  or 
absence  of  cooling,  and  the  noraallsaA  heat-transfer  rate  is  even 
less  so.  Tbls  is  true  in  spite  of  the  fhet  that  there  is  a  Sharp 
increase  in  the  slqpes  of  the  veloei'^  and  temperature  profiles  near 


tht  wU  for  the  vaxleble-prover^fluld  etm,  hut  thiere  le  a 
eorrespoadlagly  lars*  decreaee  In  viaeoaity  aad  heat  eoaduetivl'^ 
in  the  eoUlFgaa  layer  near  the  nalL.  Sie  net  effect  la  a  nodarate 
deereaae  In  ahaar  for  aaall  W  (!•••«  atroog  eoalliBg)i  aad 
an  Inalgoifleantly  aaall  deereaae  In  haat-tranafer  rate. 

liaeuaalon  of  Beaults 

Aaong  the  correction  effeeta  that  eere  ecBaldared,  the 
behavior  of  the  dlaplaeenent  tern  parallela  that  of  the  boundary^ 
layer  teas,  aa  could  be  e^qpected  fron  the  close  relatlcnship 
between  then.  Both  the  shear  and  heat-transfer  rate  are  increased; 
this  is  in  accordance  with  the  sign  coaoventlon  that  was  adopted, 
which  iagdles  that  the  Increase  is  due  to  the  increase  in  velocity 
gradient  A  (it  nay  be  zioted  that  in  nost  eases  t\  will 

decrease  due  to  the  displaceaient  effect  of  the  boundary  layer,  thus 
the  sign  of  the  dlaplaeenent  coefficient  will  be  negative,  aad  in 
reality  there  will  be  a  decrease  in  shear  and  heat-transfer  rate. 

Tor  very  strongly  cooled  coovressible  boundaxy  layers,  it  la, 
however,  possible  to  have  negative  displaceoeat  thicknesses, 
lagplylng  an  apparent  "shrinkage"  of  the  body,  thereby  increasing  /\ 
and  the  shear  and  heat- transfer  rates  likewise.) 

The  curvature  effect  tends  to  decregse  the  shear  at  the 
stagnation  point;  this  decrease  becoaws  oaller  for  swell  cooling 
ratios,  nw  heat-transfer  rate  is  affected  differently  by 
curvature  for  the  tw»»d1aens1flna1  and  axially  syaaetrlc  stagnatioop 
point  flows.  In  the  foxner  case,  the  heat  transfer  is  deereased,  aad 


•ijBOft  eoagplAtcly  unaffected  liy  tlie  eocUng  ratio.  On  tbe  other  hand, 
for  the  axially  aynnetrlc  caae  the  hBa-^txaiiefer  rate  is  IzKreaeed 
hy  curvature;  thla  Increase  heccnes  ■atller  for  aeall  cooling  ratios. 
Ihere  Is  no  slaq^e  explanation  for  these  eiuvature  effects;  they 
arise  fron  the  Inhonogeneous  tenu  la  the  eoppanded  differential 
eq]uatlans,  and  also  from  the  modified  boundazy  condltloos  in  the 
invlseld  "outer"  flow.  For  example,  it  Is  apparent  by  Inspectioo 
of  the  velocity  profiles  of  Figure  4,  that  the  negative  slope  of  the 
Invlseld  profile  has  an  effect  on  the  entire  vlse<ms  layer.  Another 
effect  Is  the  different  pressure  gradients  ejqperleneed  by  adjacent 
layers  of  fluid  in  the  viscous  layer,  due  to  the  eeatrlfugal  pressure 
rise  across  the  thick  curved  layer.  Oiere  are  aaijy  other  terns  in  the 
BavlezwStokes  eqaatl(»is  that  contribute  to  this  effect;  no  atteapt 
has  been  made  to  identify  them  separately  since  the  effects  almys 
occur  simultaneously. 

QLzeet  caaparlson  of  these  curvature  results  with  other 
theories  oaxmot  be  made  for  the  sphere,  because  the  fblly  viscous 
shoek>layer  theories  (references  13^  and  31 )  include  the  displacement, 
curvature,  and  vorticlty  effects  slmltaneously,  and  for  the  shock 
denal'ty  ratios  employed  the  vortleity  effect  predominates.  However, 
for  cylinders,  idiere  the  vorticlty  effect  Is  of  second  order, 
Hoahizaki's^^^^  theory  shows  an  increase  la  haat^traaafer  rate  at  all 
Reynolds  numbers,  even  the  large  ones  (where  presuaeUy  the  secondp 
order  effects  should  be  Insignificant).  Ibis  is  contrary  to  the 
predictions  of  the  present  analysis,  according  to  which  the  first-order 
corrections  that  were  implied  In  the  fUlly^viseou^layer  theocy. 


lUMly  the  displMaaent  and  curvature  ef facta,  are  both  negative,  i.e., 
tbejr  vUl  tend  to  reduce  the  heat  txaaafer  Inatead  of  Ineroaaing  it. 
bpanaloQ  of  Eoatalaald'a  eguatlcna  in  the  eicpaaaica  pamneter  of  the 
preaent  aoalyaia,  ^  ,  reaulted  in  the  aaae  eqjtaationa  that 

Here  derived  in  Chapter  U,  above.  It  me  not  poaaihle,  therefore, 
to  readve  thia  discrepancy  between  the  theories.  CoaqparlBoa  with 
the  eiqpexlnents  of  Tewflh  and  (tledt^^'^^^  se«Ba  to  inp3y  that  a 
redueticQ  in  heat  transfer  due  to  curvature,  rather  than  an  increase, 
is  in  better  sgreenent  with  e^erlaental  results.  (A  note  detailed 
disexiasioo  of  cosqpailson  with  experiments  foUowe  in  a  subsequent 
paragraph.) 

In  agreesmnt  with  previous  reports  (references  13,  19>  31  > 
etc.),  the  existence  of  vortlclly  in  the  invlacid  flow  (to  the  first 
order  present  in  the  axially  syanetrlc  case  only),  tends  to  increase 
both  the  shear  and  heat- transfer  rates.  33>is  increase  is  appreciably 
larger  at  snail  cooling  ratios. 

Oie  behavior  of  the  term  due  to  slip  and  tenperatiue  Junp 
at  the  wall  is  especially  interesting.  The  eag^slom  parameter  for 

A 

this  tern  (cf .  Chapter  II)  is  essentially  ,  the  mean 

free  path  at  the  wall  divided  by  the  boundary^lsyer  tlilclmess 
pazesKter.  The  term  contains  two  separate  effects;  one  effect 
arising  when  ^  other  arising  due  to 

the  difference  between  these  two  constants.  Considering  the 

K ,  '  V<  2_  effect  first  it  has  already  been  noted 

(cf.  equations  (2.39)  and  subsequent  paragraph)  that  this  term 
leaves  the  wall  heat-transfer  rate  unaffected.  Slollarly,  using  the 
locown  identities  for  the  correction  functions  (e.g.  4.4)  in  the 
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exprsstlon  for  the  aheer  correction  (^•3)>  end  oonpexlag  eUli 
houndnxsr^layer  eqjuetioQ  (2  >63)  one  lonedlntely  oibeervee  that 
at  the  wall  Ix^epective  of  VJ 

The  other  slip  and  tsnperature>JuBqp  effect,  the  texm 
proportional  to  (i'^-  ,  affects  both  the  bhear  and 

the  heat-tzaasfer  rates  at  the  waJl.  She  shear  is  inersased;  the 
increase  varies  from  0  to  large  values  as  W  varies  from  1  to 
sonll  values.  (The  uzicooled  nail  also  represents  the  constant-property^ 
fltiid  ease,  in  this  case  can  have  oibvlaualy  no  effect  on 

the  shear  since  the  mcettntua  and  energy  eqpations  are  not  coupled). 

But,  as  V  becoBies  saadl,  the  osan  free  path  near  the  wall, 
iMch  determines  the  amount  of  slip,  and  iddeh  appears  in  the 
expansion  pstrsmeter,  also  becosMs  ssmll,  since  the  cooled  gas  near 
the  wall  is  more  dense.  It  is  of  interest  to  find  the  cosdilned 


effect  of  strong  cooling  aa  the  two  ccnoetlng  effects:  the  increased 


correction  function  and  the  decreased  naan  free  path.  Since 


one  may  write 


hi  yy  —  l>  -L"  'fc  ^  h  hy  ^  ^  ^ 

^  f.  ki  &  Ai  IT./  ^  ^ 


The  quantity  is  plotted  in  Figure  6; 

the  graph  shows  that  as  W  changes  from  1 .0  (no  cooling)  to  0 


(strong  coding),  the  qjoantity  increases  from  0  and  then  appears  to 

a 

approach  a  constant  vdue  for  very  swan  W  The  coniblned 


effect  on  shear  of  the  "slip  and  Juop”  texm  is  then  an  increase  for 
strong  cooling,  because  as  VI  0  ,  the  'ty  tern 

dominates;  Whereas  in  the  region  near  W  I  the  only 

slgDlfloant  effect  is  the  tern  -  I  ,  and  the  shear 

is  thus  decreased. 


The  MLll>hflat>trttasfe]>zmte  eorreetlon  efftet, 


sb0M  •  b^vlor  BlMiT  to  (0)  Ow  effect  tends 

to  decrease  tbe  heat-txensfer  rate;  as  W  ’becosws  very  saall 

the  decrease  becones  excessively  large.  But  if  the  qtaantl'^ 
is  plotted,  it  resislas  reascoahle  in  magnitude  throucd^t  the 
entire  range  of  W  Ihese  observatlOBS  indicate  that 

for  the  affect  that  is  under  caoslderatlon,  namely,  the  correction 
ailsing  due  to  the  difference  in  the  eonstaats  of  proportionality 
for  the  velocity  slip  and  temperature  Jump,  K,  and  ,  the 

proper  parasKter  that  deteaines  the  order  of  magnltiide  of  this  effect 
is  not  'y'l  $  as  uas  formerly  supposed  (e.g.  references  9t  31  > 
etc.)  but  rather  Ihls  liqplies  that  the  reduction  in 

heat  transfer  at  the  stagnation  point  of  a  blunt  body  due  to  slip  could 
be  significant  even  for  the  osise  of  strong  cooling  in  hypersonic  flow. 
JUst  how  large  this  reduction  any  be  in  different  flight  regimes  is 
indicated  by  the  lines  of  constant  plotted  in  Figure  3« 

Conpsurisan  of  Besults  with  Biqpezlaents 

It  is  finally  of  Interest  to  coopars  the  nowrlcal  results 

presented  in  this  chapter  to  experlamtally  deteziained  properties  of 

loi^Iteynolds-number  stagnatlan-polnt  flow.  For  a  partic\>lar 

ejgMrIneQt,  It  is  accessary  that  the  appropriate  flow  psunaeters, 

Bsynolds  number,  stagnation-point  velool^  gradlant,  mean  free  path, 

etc.,  he  knewn,  and  then  the  resib.ts  of  tbe  present  analyals  oan  be 

applied  to  predict  the  flow  properties.  In  one  series  of  er^exlsients 

(26) 

NSloe,  Butowskl,  and  Chsn' 


measured  beat-transfer  rates  at  the 
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blunt  nose  of  a  heBlBpbere  cylinder  placed  Into  a  low^denalty 
byperaonlc  shock  tunnel.  But  these  expezlnental  reeults  are  not  too 
well  suited  for  coigparlson  with  the  present  analyele.  Ibe  reason  for 
this  Is  partly  the  difference  between  the  fluid  properties  of  the 
present  analysis  and  those  of  the  high  tenperature  dissociated  air 
of  the  shock  tiinnel.  Tbls  difficulty  could  be  overcone  Igr  taking  ratios 
of  the  theoretical  loi^iteynolds>number  and  boundazy* layer  results;  and 
then  applying  this  ratio  to  available  dlssoclated-air  boundagty- layer 
solutions  to  obtain  a  reasonshle  theoretical  prediction;  which  could 
then  be  cosgiared  to  the  e:iqperiaental  results,  (ihis  was  the  scheme  used 
by  the  authors  of  the  e^g^rimentS;  who  cosqpsuwd  thdr  results  with  the 
constan-Ufluld^property;  "exact,"  viscous  shock-layer  solution  of 
referanee  13«)  A  more  serious  difficulty  in  using  these  experiments 
as  a  basis  of  eonpaxlson  is  the  predominance  of  a  veiy  large  vortlcity 
effect;  which  puts  the  results  beycs^  the  reasonable  validity  of  the 
expansion  procedure  of  the  present  analysis. 

For  cylinders,  an  extensive  series  of  low>neynolds-nuid9er 

(JjO)  (41) 

flow  measurements  were  performed  by  Tewflk  and  Chedt'  .  Oiese 

experiments  were  performed  in  the  ltocl>*nuBd>er  range  of  1 .3  to  $.7> 
and  with  low  tenperature  (and  therefore  non-dlssoelated)  air.  Ohus 
the  fluid  properties  that  were  used  in  the  present  analysle  are 
exactly  those  that  sure  applicable  to  these  experimental  conditions. 
Furthermore;  for  the  range  of  flow  paraneters  that  were  cnployed 
in  these  tests,  all  sigalfloant  Icn^Bsynolds-nukber  effects  are  about 
equally  large,  sufficiently  large  to  be  laportant,  yet  not  large 
enou#  to  preclude  reasonable  applicability  of  the  ejqpanalon 
procedure  of  the  present  analysis.  For  these  reasons,  a  eonpaxlson  of 
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the  eaqperimenta],  results  vlth  the  present  theory  seeaed  espeolsUy 
appropriate,  and  was  undertaken  in  detail.  Ibe  pertinent  calculation 
procedure  Is  described  In  i^ppendix  E;  the  results  of  the  calou^tlons 
are  presented  In  Teble  IV. 

The  free>streaa  Nsch  nuniber  M  oo  >  Bsynolds  zumber 
cO  (ranging  from  37  to  4l00),  and  the  waU-t>*etagaatloiv-teaperature 
ratio  Vj  (raaglsg  from  0.24  to  0.74)  are  the  Independent  parsmeters 
ehlch  are  the  Input  necessary  for  application  of  the  theory.  In 

addition,  the  stagnationrpolnt  velocity  gradients  were  also  measured 

A 

and  recorded;  from  these  the  slgilflcant  Beynolds  nraher, 

(based  on  stagnation  fluid  properties)  could  be  calculated.  Ihe 


Isnrerse  square  root  of  this  quantity  is  the  er^ansion  parameter  of 
the  present  theory;  it  detemines  the  size  of  the  curvature  correction. 
Both  of  these  quantities  are  tabulated  In  the  table,  the  former 
ranging  from  25  to  990>  the  latter  from  to  20^.  The  quantity 
that  is  significant  for  the  velocHywaiip  and  temperature- Junp 


effect, 

1  .Qi  to  14)1. 


is  also  tabulated;  it  ranges  from 


Baaed  on  the  above  infonaatlon,  the  theoretieal  results 
discussed  earlier  in  this  chapter  vere  used  to  predict  stagnatioor 
point  heat-transfer  rates.  Predictions  were  based  both  on  the 
constant- flvcLdr property  and  the  varlable-fluldpproperby  theories. 
Curvature  and  velocity>8llp-teeqperatura-JuBp  corrections  were 
considerad.  Ho  other  low'Bsynolds-nusiber  effects  oeeur  sinoe  there  is 
no  first-order  vorticl'^  correction  for  cylinders,  and  the  effect  of 
boundary>layer  displacement  on  the  external  flow  has  already  been 
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accounted  for  ty  ualag  the  expexlneBteUiy  aNMorod  reloeitgr  gnftleate. 
IBie  predleticae  ere  presented  In  terns  of  a  eosoarlson  to  tbs  hea^ 
transfer  rate  based  on  uneozxeeted  eoastaat>fluld*propertgr  boundsuy^ 
layer  theory^  viz. 


(D  . 

^Ou%c. 


0  s\th 


(\  -  w)  l-vTiy^ 

All  the  eorreetioas  that  were  oaleulated,  are  shown  in 
mhle  IV>  as  a  fraction  of  the  shove  nuaher>  In  the  foUowing 
order:  inccapreealhle  curvature  correction^  IneoavreBsHile  taaperatuxe> 
Jusp  correction^  total  Incoevresslble  correctlon>  conpressihle 
boundaxy^lsyer  correction>  coz^resalhle  curvature,  tengperature-Juqp, 
and  total  cosQiresslhle  corrections.  Hie  eiqpezliiental  results  are 
presented  on  the  saDe  basis,  as  a  deviation  from  the  prediction  of 
Incoavresslble  boundary-layer  theory. 

Hie  tabulated  predictions  Indicate  a  redaction  in  heat 
transfer  due  to  aU  the  effects  that  were  considered;  the  reduction 
Is  lars^r  In  the  ccnpressible  case,  especially  for  the  teoperature- 
Juap  term.  Eowever,  the  experimentally  naaaured  heat-transfer  rates 


uniformly  differ  from  the  Ineosq^resslble  boundary-layer  prediction  by 
significantly  larger  asiounts  than  either  eosdlnatlon  of  correction 
terms  predict.  In  order  to  investigate  the  posslblll'ty  of  a 
relation  between  this  discrepancy  and  the  dlsplaeoient  effect, 
displacement  coefficients  Cp  were  calculated  for  all  the 
points.  Hie  calculations  were  baaed  on  comparing  the  experimentally 


obtained  velocity  gradients  with  the  hldhrBqniolde-aRndier  velocity 
gradients  given  by  Rsahotko  and  Beckslth^^^  (of.  Appendix  H.)  Hiere 
appears  to  be  no  relation,  nor  do  the  observed  diserepanoles  show  any 
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otter  olavloufl  regularly.  Hat  dlserapnaiM  tten  rwin  uauq^ljdiMd, 
pending  further  eoaperlaoa  with  other  eaqyerlaentel  zemlta.  Bsverthe- 
leu,  It  oaa  at  least  he  ascertained  that  the  pzedUtione  of  the  theozy 
Bhav  the  same  tread  as  the  experlaaotBl  results;  anaely  a  rsduetlaa 
in  heat  tzensfer.  This  is  slfpiifiaanct  titeo  eoepared  to  the  prsdietlons 
of  references  1 3  and  31  >  which  show  an  increase  in  stagnatl(ne*pclnt 
heat^transfer  rate  for  a  cylinder  at  lav  Itaynolde  naaters  (tteas  ealcu- 
IcatioDS  wealweted  the  decreaae  due  to  velocity  slip  oad  toeq^rature 

JUBV.) 


APPmX  A 


tezlTmtloa  of  the  Squatloos  of  Motion 

The  Ibvlei^Stokea  (aonentun)  •g^atlona  for  atendy  flow  eaa 
be  vritton  dom  In  vector  foxn 

g  ^  ^  i  ^  s.  ^  ^  J>b«r  ^ )  4- 


4.  X  -■  *^1 

A  general  orthogonal  coordinate  ayatea  can  be  defined  by 

the  dlreetlona 
^ 

with  eoordlnatea 


and  natxle  funotlona 


(A.2) 


^  At. )  s 

Let  /,  =  conatant  define  the  planea  of  ayaaetxy  of 
the  flow  field)  ao  that  the  deilvAtivea  ~  of  all 

qMaatltlea  vanlah.  For  thla  caae,  a  eonpreaalble  atraan  funetion 
can  be  defined  by 


r  A  .  0  Hi 


(A.  3) 
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iMoh  IdntlMlly  satisfies  the  oootliBiltgr  e^mtiai  for  steadgr  flcnr. 
Substltutiag  (A.2)  saA  (A.3)  into  (A.i)  sad  uslag  veetor  elgehra^  the 


e. ,  eoB^oneat  of  the  aossatus  equation  eaa  he  vrlttea  dona 

as 

_i— \  .  i p  vv»5v  .-i  h. '3«»  \ 

^>1  **•  *’'•  ‘‘i  3’'>'  ”  vJ,  ■*'  SbQEii  Tp^  )  ~ 

_  u  u  "l!  f/ j.  'Jj? f  4./ j-  ’jf  f  1  ^  ,  r  '*j'  1  -1-  ’jf  .  v»,  St  j_  ^jr  1  . 

4-  U  i  -  1  f  1  + 1  1  -i-r  - 

>»ij  3  L  3»J  3.,  i"!  ».^j  K,  3»«,  iVAv^iWl  9X> 

-"-S'  ^  1  _! _  \1  K.kv  3>^,  .  1  ViAv  tbl*. 

Hj  L’n  "k:  >*>  S  >'‘>i 

4-il^-i  ^  2  Ha  0<r  1  .  I  ^  ^  ^  1  - 

i-  I  k  (  2.  ±±t  ^3*'  +  2  ilk  1  /»  J.  X 

ainllarly,  the  coapcaeat  heeosMS 

-  >  -  /i  f'nt  '^J  2i  +h‘  ii  ).  i  fl  hi,  + 

\.''v^'5  ''1  ’'**  S  Vj  5xi\^»i  k,v*i  a*i 


1  k  \  L.'lsr/'i  kr\^jL 
3xjV,^V)j_  5-jJ  3xjLv^>x  3x,  J  VWki^xjM 

r  3_±  X  —  ®-*  +  b  '2-  "1  ^  Vv  k  ^vi¬ 
sta-,  ax^  V,Vv^,  v»j  axj  3«j  Vvb ’>■3  \  ^  '  *■  ■a'»3  3 

.  4  i  2  r  Xt  \i  - Xi  n  M  i  _J_ 

■  3  Vij3Xi  s’'Vv,Vik3  L^-'J  ^’'5-J  Vj3xj  S^'kb' ■=»>'.  »-j, 

-  li  3^  \  _  -!-  XT  — ! — r  \  }■  kb  +  2.  bk  1  _ 
3"!  '  V>,l,t  s''«'"vki  L^x,  \,,  >>',  >Xj  Vij 

'  3j-r^  Vv.  ^  3  WJ  3J|>1  _  -1—^  U:  PJ:  U  ■ 

'  Vi^Wj  L^^SKvki  axj  'ik.Vvb  I  “1  "J 

Sit  ii  k  X  3L  k  fi  kt  ^  k  >-b  I 

3’">  »«i  3’'>  *.bV\x  t>x,  J 


(A.  5) 
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Vor  tha  ■pedal  ease  of  a  "poOar"  oootrftlaata  ^toe  tlie 
coordinates  end  aetxle  functions  are 

x,*r  j  x^= 

(A.6) 

vbere  ^  0  Is  the  plane  tir(>>dljaensloaal>  and  n  I  the 

arlally  syaaetzlc  ccum.  (Qtese  are  the  connreatlcnal  cylindrical 
and  spherical  polar  coordinate  aysteas  respactlTely) .  Using  (A.6) 
in  (A.4)  and  (A.$)>  end  perforalng  all  the  algehrai  the  two  acsMntum 
equations  hecosM 

_ _ ! _  I  yr^»  _  "il  _  Y\  Sy- 

5  fr 5  V-  yV  ^ 


+ 


y»-Vr»  .  o^V'-)y»  7^0- 
r"*-  i 


» 

(^r 


I  -  A-  - 


V  Ci->v')Mf»^Sr  ^  1  1  y;>Sr8»-  ^  T+vj  '♦vA_S^  2L 

i  S'*"  i  5  5  i  S'^  'f 


- T —  — ^ ^ 

J.  ^  ,  xAOr>&»W>«.  _  ^ 

.  ^  y>«-  " 

y  ^  y'-  S 


t»-  ^  Vyy-  a*-  ^  'Yfk-  ^  It  .  ^Y /*»  St»  ) 

V  %  S  J  J 


(A.7) 


la  ttae  'f'  dlzaetioa,  uA 


,  P?  _  _  /i  y*^h*iT  i£'K.>^*5» 

■»'  r  jc«rfw.sri'^  "  3  5  •‘- 

H-  I  +  t  ^  ^  i  ^  ^  - 

-2  y-  .  ft  _  '<'«.>.<&►  a.  Vv  Y»>*-y 

i  •A>  4  »*•  4  s’-  '>i-w »• 


K«r  ^  2.  Hr,»  _  H[»»»  ^ 

r"’-  «■  ”  *’ 


_  Y  >*-»  !«•  _  2-»v»  >>•«>  4.  ’».»»  hi  _  2- _  V\  4. 

S  y'  fr- 


'r'^nA-.O'  '  '  • 


ST' 


-  Cv-V'"! 


_  0_iJ2}2Ci2il  -  Yt  ^  2  Yrh-Sf^  1 

;  T-  "  r*-  s  3?-  J 


(A.8) 


in  the  ^  dinctloa.  33w  mbaeripts  v  aaft  d-  danote  the 
resp*etlv«  partial  dexlvatlTes. 

lov  the  energy  equation  for  steady  flov  eaa  he  eonsldered. 
In  vector  notation 

3  .  (^^«-(v/.T  -  :SL.  t^-iy  Z'-  -Y 
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ttqpmdlag  in  teau  of  tbo  gnoxal  eoooriiiMtta  qrst«i  of 
(A.2),  the  oaezgy  aqiaKtloa  (A.9)  beeones 

^  r  U 1  — i _  pj:  if  1  _ 

3-;  J  La>^,  i>K4 


-1-  [1. 

■  v.v.^J»« 

jT  "i  Vit^ti  0  ^  “ 
i«,  »>«>  V»3  ^^-1, 

!*>-■ 

3. 

,  V,  dx, 

.  %.  \\i-  2- 
3Kt,  V-i 

(i 

■*■  Cw.vi)  [>\  ^ 

jii- 

|V^-^, 

»'l  \ 

.  i 

r  2  ^ 

\.3»,  sW,  a)*! 

U.V.,  ( 

r  0  'Vi 

5kV, 

h  ^ 

J2i. 

w 

1  b.  -  f  3. 

V, 

'iJt  ^  1 

2Jr^l 

^  ax, 

sV.Vt 

iXj  ><c, 

S'*>v 

2. 

[VL  i. 

_ 

1  1  JL-  1 

'ix  - 1  - 

•  I  » 

i 

1 

L»*iS 

ax,  ?x,  ! 

z 

~  3 

cw.v,. 

.  ^’'j.  5-, 

-  ■ 

^  D*>l 

i‘i 

(A.  10) 

B»r«»KHng  in  the  "poOLer”  eoordlafttes  of  (A.6}>  the  efuetioa 

beeones 

S  L%\  '  'V  JO  "  i  Pr  -  ^  Ve- )  =  + 


^  ^  O-'  ^  r"-'  \ij»  ^  ^ +  ^  J 

(A.n) 

lAere  ^  le  the  dlailpetlan  funetloa. 
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for  future  referenee,  the  expreeeloB  for  la  the 

generel  eoordlmte  eyeten  le 

-  :  r  '«>  -3^  \8  ^  u,  ^  ^  ]  —  t-  11  4.  ^  ^ 

(A.  12) 

for  the  special  ease  of  the  "poLar"  coordinate  qrstem  of 

(A.6) 

il  -  ^ 


^  v\  UftV  ;  »»»■  _ 

'T^'ki.w  S’ 


(A.  13) 


ACTBQgl  B 


terlvatlon  of  Boundary  Condltiooa  at  ttai  Sodldpaas  Interfile  • 
Tbe  Boundary  eoodltloBa  for  a  gu  flcMli^  along  a  solid 
body  ean  be  derived  from  the  kinetic  theory  of  gsses^  and  are  given 
by  osny  authors,  e.g.  Schaaf  sad  Ghnoibre^^^^as 


0(0')  i 


'au. 


1-0 


.  i  K  ^ 

(V  sT  3,f 


pi*.  7+'  %.  ''4*-0 


(B.1) 

(B.2) 


ehere  the  sane  reference  glvestthe  mean  free  path  la  terns  of  the 
fluid  properties,  as  foUom: 


)  ~  fz  ^  _ r  h- 

^  ^ou  ""  hsi  fjf 


(B.3) 


In  addition,  the  no  through*  flow  boundary  condition  can 


be  used 


VJ-  to)  =S  0 

(BA) 

Iheae  expressioas  ean  be  applied  to  the  stagnatloopolnt 
flew  that  Is  being  considered  by  subatitutiag  into  them  changes  of 
variables  and  ejqpaaslons  (2.1)  through  (2.7).  To  find  0  (2.2)  and 

(2.4)  can  be  substituted  Into  (1.4). 


(B.5) 
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fiL£fttx«atl«tlag  OHM  on*  bbtalai 


ris  + 

3'^  R  1  K  L  K  H  J  K 

1  -1 1*  K  Vi  y  2.  -V  'V*  *V>  4.  . , ,  1 

+  -ir  ] 


(B.6) 


8iJBllarly«  tha  expx«iiion  for  o'  taeonss 


,  _ _ ff*  .iF5if^'  1 

- 1^7 


(B.7) 


tt. 

T)« 


tern  In  the  Blip  veloel'ty  1b  of  order  , 

and  benee,  vUl  not  eontzihute  to  the  order  1  and  older 
that  are  being  cooBldered.  Bqpatlon  (2.3)  can  ncer  be  need  to  find 


•)T 


0  I  _  T\ 


(B.8) 


and  flnnlly  egiuatlooB  (2.3)1  (2.U)  and  (2.7)  to  expreBs 

1  ( !i:;  ■*  'r^  ^  1  f  ^  T  — 

1  K,  K  V  M  L  V, 

_  kof  V.  .  '  I 

v.r)J  J 


(B.9) 


UBlng  (B.5),  (B.6)  In  (B.1}»  and  eqpatlng  like  poerere  of 


1  ly 

Heynolda-raimber  paznaeter  ^ 


ooB  obtaluB 
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.0 

Wol 

aad  (uBlag  this  result) 

i;<„  4^  nf  ^ 

auollaray,  using  (B.7)  In  (B.4) 

\tLo)  =.  0 

=  0 

W^*Uor>  (B.8)  end  (b.9)  exe  used  In  (B.2) 


(B.10) 


(B.n) 


(B.12) 


»  Is'  £  W 

^,,0)  =  ^-2^*-  [t  ^  iJt  1  ^  , 

'12  rCi; 

aince  (B.13)  sbom  that  Ao^o)  =  A- 
expresslan  oen  he  renitten  as 


(B.13) 
,  the  last 


K,o, 

•<..  Je  )€?;  I  \  c,^  (w 


(B.14) 


APHWP  C 

■unzloal  Solutloi  ot  lovlMid  Qypanonle  Hair  iMur  8t«0»tlon 
Points  of  Siptaare  and  Qylindtr. 

OSie  first  step  In  the  miaerLeal  sOlutloB  Is  the  solution  of 
eqjuntlon  (3*7)«  staleh  relates  ^  to  e  .  For  the  ease  of 
the  eyllnder«  the  solution  ean  he  ohtalned  graphioallF'  ‘by  wilting  the 
eqiuation  In  the  font 


K, 

T 

1 


The  functions  on  the  two  sides  of  (C.1)  ware  found 


(C.l) 


nueerlesUy  as  functions  of  their  xwqpaetiva  axgianents  with  the  aid 
of  Bessel^funntlon  tahles  and  then  plotted  on  Ike  ssbb  gmph. 


Oorreopcading  values  of  the  two  aqgusmts  ware  thsei  found  gmiMcdlly, 
fm  which  the  corresponding  values  of  —  and  e  ware 
eoaputed.  These  results  ware  then  used  in  (3.10a)  to  oaloulate  the 


veloeity^gradient  paxaneter. 

For  the  sphere  a  slinllar  poroeadure  was  folloiiad>  except 
here  a  graphical  solvtlon  was  not  nasesoaxy  sdnse  (3.7)  is  a  qiuadratlc 
in  6  >  and  can  he  solved  aaalgrtiaally 


6  -  do( 


e  = 


<100 (it if t 


a  vF 


(C.2) 


The  result  of  this  sdution  was  again  uaad  in  (3«i0h)  to 

chtodth  ^  and  V  R  for  the  sphere. 

■u 
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APPHnnx  D 


CkleulatlooB  of  Flour  Flameters  at  Varloua  Fllc^t  ElpeedB 
and  Altltudea  In  the  Ataoqpliere. 

In  order  to  calculate  the  flow  paraaetera  that  are  of  Intereat, 
the  following  three  quantltiea  have  to  be  calculated  at  varloua  fillet 
apeeda  and  altltudea;  klnenatle  vlacoalty  and  aean  free 

path  A  5  ,  both  at  the  (lanrladd)  atagnatlon  point,  and 

(for  atveraoole  fUc^t)  e  >  the  bhoek  denalty  ratio.  Itae  altltudea 
were  choaen  at  50,000  ft.  Intervale,  atartlng  at  1^,000  tq?  to  and 
Including  400,000  ft.  The  fU^d^t  apeeda  were  grouped  Into  three 
reglnea:  .C^)  perfect-gaa  regime,  10,  100,  1,000  ft/aee.,  (2)  non> 
dlaaoclated  regime,  3*000  ft/aee.,  and  (3)  real-gaa  regime,  7,000, 

10,000,  20,000  ft/aee.  Oae  ABDC  aodel  atmoaphere^^^  waa  uaed  to  find 
the  freo>atreaa  denalty,  teaperature,  and  q^eed  of  aound  (the  latter 
by  extrapolatloa  at  the  two  hl^ieat  altltudea).  The  atandard  reference 
temperaturea  and  denaltlea  of  7~  atandard  -  3l8.69*B  and 

^atandard  -  2.3769  x  10  ware  uaed. 

Different  calculation  procedurea  were  uaed  In  each  of  the  three 
reglnea,  aa  foUowa. 

In  the  perfect-gaa  regimB,  the  frea-atraam  Mach  nudber  waa 
detezmliied  flrat;  then,  ualng  Y:  1.4-  ,  atandard  compreaalble- 

flow  tablea  (e.g.  reference  l)  were  uaed  to  find  atagnatlon  thermo¬ 
dynamic  propertlea. 


,86 


87 


At  3iOOO  ft/sec.  KauAuan's^^^^  tables  (based  on  the  Beattle- 
Brldgenan  eq]aatloa)  were  used  to  detexnine  the  quantities  of  Interest. 

The  tables  give  shock  density,  tesipenture,  and  pressure  ratios,  shock 
Mash  number,  stagnation  tenperatiire,  and  press\xre  directly  as  functions 
of  free- stream  Mach  number  and  altitude.  Perfect-gas  relations  were 
then  uded  behind  the  shock  to  detenglne  the  stagnation  density  from  the 
known  shock  theimodynamic  properties  cuid  the  stagnatlcm  pressure  and 
temwrature. 

(12) 

In  the  real-gas  regime,  Eochatlm's  ’  chart  was  used  to 

find  the  teiiQ>erature  behind  the  shock,  and  the  density  ratio  across, 

and  henee  density  behind,  the  shock  at  the  various  flight  conditions. 

(7) 

Feldman's  Molller  diagram  for  air  was  used  to  detendne  the  effect 
on  these  thermodynamic  quantities  of  the  Isentroplc  compression  behind 
the  shock.  The  enthalpy  change  during  this  compression  could  be  determined 
since  the  velocity  after  the  shock  was  known  from  continuity  considerations. 
(Some  care  had  to  be  exercised  in  using  this  multiplicity  of  charts, 
since  both  Hochstlm's  and  Feldman's  reference  conditions  are  slightly 
different  from  those  of  the  AIQX!  atmosphere).  Compressibility  factor 
z  at  the  stagnatlcxi  point  could  also  be  determined  from  the 
Molller  diagram. 

(10) 

Hansen's  calculations  were  used  to  determine  the 

viscosity  at  the  stagnation  point  as  a  function  of  tenperature  and 

(at  high  temperatures)  of  pressure.  (A  graph  of  stagnation  pressvires 

at  various  flight  velocities  (uad  altitudes  Is  given  In  the  same 
(37)  (25) 

reference).  Schaaf  ,  and  also  Maslen'  '  use  the  following 
expression  to  calculate  the  mean  free  path 
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«h0r«  ^  ud  T  ar*  la  units  of  ft^/soe.  aad  *R  rospoetlTely. 

Constsat  1 .255  =  sboro  Is  Issa  than  2^  diffsrsnt 

froB  1.278  «  — ^  ,  iMeh  Is  used  tgr  Hsasen^^^^  aai 

.  5  fli? 

Pattersoa^^^ . 

USlag  the  shore  proeeduzes  to  detezalne  , 

sad  ^  ;  and  using  Figures  1  and  2,  all  the  pazaasters  given 

la  eqjuatlan  (3*1^)  could  now  be  calculated  at  the  various  altitudes 
and  fU^t  speeds.  Ihe  curves  of  Figure  3  were  then  tibtalaed  bgr 


of  cross  plots. 


APMnBX  a 


Calculation  of  PoMr  Lava  for  VhrlaUa  Huid  Frqpartleai 
(10) 

Hansen 's'  '  tables  vers  \uied  to  detemlna  the  specific 

heat>  viscosity,  and  heat  ccoductivity  aa  functions  of  teaqperature. 
Stese  results  vere  plotted  on  log-log  graph  paper,  and  tangents  to 
the  curves  vere  dravn  at  various  tenperatures.  The  results  vere 
aa  follovs: 

For  specific  beat,  <<  ' 

t<  mO.Mik 

For  viscosi'^,  <  1 

at  5001  *  0.661 

at  10001c  00  =  0.606 

at  aoooiC  oOsO.^ 

For  beat  conductivity,  Ai,  e.<  T 
at  500*«  t-0.715 

at  lOOOlC  ^=0.702 

at  15001c  e=0.67^^ 

The  Frandtl  nuaber  at  four  representative  toaperatures  is 

given  aa 


Tic 

Pr 

500 

0.738 

1000 

0.756 

1500 

0.767 

2000 

0.773 
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OoMlteriag  1200*K  u  tlw  dtilz«bl*  aidlaa  t«p«zmtur«  for 
tb«  nago  under  eoneldAration*  and  aaglaetlag  tha  hmU  variation  in 
Prandtl  BUHibar«  ttaa  following  propartiaa  raprasant  an  optinun  "fit;" 
Pr  ■  0.76 

ei  ■  0<11>  £>  0i69 
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utmax  j 

OtrlvatloB  of  EqiuatiOBUi  for  Constant  Fluid  Frppertlas 
lot  tha  ease  of  constant  fluid  propartiaa,  the  analyala  of 
CSuvtar  I  aad  II  Is  modified.  !B»  perfect>gas  eqiuatlon  of  state  is 
replaeed  'by  3  =  constant >  aad  the  temperature-dependant 

Ylscosity,  heat  conductivity,  and  specific  heat  are  replaced  by 
constants.  Thus  there  vill  be  no  eou]^ng  beteeea  the  naaentuB 
and  eaeny  eqioatlons;  the  acawntuai  aQoation  can  then  be  solved  first, 
aad,  using  this  solution,  the  energy  equation  subsequently. 

In  the  ixnriscld-floir  solution  of  Chapter  I  result  (I.9) 
vas  derived  from  the  energy  eqi^ation  only,  aad  hence  resMlns 
unchanged;  likewise  results  (1.13)  thxou^  (1>2C>)  sn  all  derived 
from  the  mcaentuu  eqiiations  only,  and  therefore  also  remain 
unchanged.  The  former  detemines  the  first  term  in  the 
temperature  expansion,  the  latter  aU  the  terms  that  were  considered 
in  the  streasH  function  aad  pressure  ej^anslons.  These  results  can 
then  he  used  unchanged,  as  they  are  presented  in  equations  (1 .22), 
(1.24),  aad  (I.23);  these  furnish  all  the  'boundary  conditions  that 
are  necessary  for  solution  of  the  viscous  flew 

^  ~  S  ^  +  nv)  I  +  '•«  ] 

T  ^  -  I  +  Cum  ^ 

P=  Pj-  sA' I  1-^'-  ^  J 
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Dor  tk«  vlseoua  flow  tlM  dlffiamtlal  cfurtlaM  hinr*  to  bo 
■odlflod  iB  MoordcBoo  vlth  tho  eoBotaaejr  of  fluid  proportioo 


0 

u 

)  V.  ^  ^ 

—  =1  , 

V*  =  1 

K  ' 

K  ' 

TboM  Tooults  cun  bo  uflod  throughout  Ite  aaiOyolB  of  Ghoptor. 

U}  e.g.  (2.9)  throuc^  (2<12)  boeono 

==  (^*3) 

Sw  first  three  of  these  eq^stloiis  eon  bo  Intogrotod  laMdlotoly 
M  la  (2<27)>  slBco  the  bouadaxy  eoadltions  given  bgr  (F.l)  rwsln 
unehoaged 

=0 

|ji  '  (P.4) 

Using  (F>2)  sad  (F.4),  bouadary^logror  oqioatloas  (2.14)  oad 
(2.17)  boeoae 

,  1  I  '*  1 '  *•  1  I  '<1  ^ 

C<  +  »')'^oA:o  -  >T«  +  I  +  Y®  “  ® 


(F.5) 


^C»  +  »»)  — -0 
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9m  bouadazy  eoodltlons  at  tiM  mil  mr»  dbtalMd  from  kloatle 
tlMozyi  and  xmnin  wunhangad,  m  cl'*^  (2.20).  fllae*  tte  loflseld 

flair,  giran  by  (V.l)«  Is  unehangsd,  (2.25)  1«  squally  ivplleabls. 
Suaazlxlag  the  boundary  eoodltlons  then 

® 

0 

^.(0)=  W 

11  1  (r*6) 

AJ  OO 

-►I 

Zhe  solutloa  of  (7.5)  subjeet  to  (7.6)  eaa  be  perfozmd 

la  tvo  steps;  the  aoMtatua  solutloa  Is  SLssMas’s  and  icaaaa's 

elassleal  result  for  the  two>dlaaasloaal  aad  axially'  syaaetzle 

(38) 

oases  respectively,  as  given  for  exaaple  by  SchUehtlag'  ' .  9m 
energy  equation  Is  Uaear,  and  caa  be  aoznallsed  with  respect  to  the 
azbltrazy  tanperature  ratio,  by  settiag 


'(  A'I)  =  W  +  (I  -  w) 

(7.7) 

so  that 

=-0 


/>l  — >  OO 


^•(<1)  I 


(1.8) 


Bm  loltttlaas  of  (r.8)  eaa  bo  mlttoa  doni  la  tho  fom  of 
tbo  foUoolag  Imtogzol 


=■ 


I 


It 


(F-9) 


TObulotod  voluoo  for 


•ro  glTon  for  •  xvago 


of  Proadtl  aubboro  bgr  Ckld«tol]i  and  fbr  the  tw>dlJMiioloBol 


•ad  oadolly  ejoaotrle  oooeo  roopoetlToly. 

She  ogMotioao  for  the  eorroetloa  texae  «•  Obtained  by 
uelag  (7.2)  la  (2.15)  ead  (2.l8).  Coaelderlag  aoMatua  oguetloo  (2.15) 
flret 


+  t  I*  as. '1^  1^(2. n 4- 1)^0 '•'a)  ^0^*  — 


I  ■t'  t  ^ 


Ueiag  (2.b7)  ead  {2.k9),  the  abore  e^^reeeloa  becoaMS 


s.  -  cw-vi')^  ]  ■'*  ioj) 


ebloh  deflBee  opexatore  ead 


(7.10) 
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Hw  rwBzk  Mdft  abonre  eoneemlag  th«  tooundaxy  ooodltloas 
for  tb*  bouDdaxy^lagror  tens  is  aqviaSJgr  appUeable  to  tbs  flrst-ordar 
eorreetlon  texaa  and  ;  bot  the  mil  and  "outar" 

boundary  oondltiona  are  unchanged  froa  the  variable-fluldpproperty 
eaae.  At  the  mil,  from  (2.35) 

^  ( O)  s  0 

Aa  for  the  ▼arihble-flul^property  eaae,  the  linearity  of 
the  eqiuationa  can  be  uaed  to  separate  the  various  "effects,"  lagged 
by  the  arbitrary  constants  appearing  In  the  boundary  conditions.  Ualng 
the  results  of  Chapter  II  then,  the  eqiuatlon  and  boundary  conditions  for 
the  curvature  term  are 

io)  •  0 
f!c  l^)  ^  0 

7—"^  (m) — 

Vor  the  slip  tern  '  I 


(P.12) 
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(W) 

{..Jo)  -0 

I  '  » " 

4-'^'-C’i) -* 0 

Bm  solutlaa  at  (7.13)  !■ 


Om  dltplaeannit  effect  le  detexaliied  "by 

-  a, 

^ 

^'5  (<>'>  =  ^ 


(F-13) 


(7.14) 


1 


— »  «>o 


ubieh  bM  the  sdutloa 


(7.15) 


(7.16) 


flaelly,  for  the  rortleity  effect  one  can  mite 
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njMUor>  ^or  thft  Toortleitsr  effcet  om  eaa  vxlte 

(o')  =  0 

4.V  (0)  ^ 

7'»~  M-'T-’Mvf') 

Iha  eorrectioa  texm  for  the  streea  funetion  eea  then  be  written 

I  4  b*  iw  y\\,  + 

'  iicT^  hi 

(P.l8) 

Ihe  eqtintlon  for  the  teagpereture- function  correction  tern  is 
obtained  by  \>slng  (F<2)  In  (2.18) 

where  the  function  ,  as  given  by(p.l8)f  is  new  a  known  function 

Obtained  fron  the  solution  of  the  monentun  equation.  Defining 


w)^,  =  .\4, 


(M9) 


and  using  (F.7),  the  following  differential  equation  is  obtained  for 
^1  ,  (pemittlag  solution  for  all  values  of  \A/  ) 


^-t-  CV4W'J  t  ^0  ■t-'V^o 


(F.ao) 
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iMoh  Otflaiti  ogpcrator 

I  u*  oWainad  ttm  (7.11) 


Bw  bomtexy  eoadltl«v  fior 


/  (7.21) 

Due  to  the  Unaoxlty  of  (7.20),  the  four  "effeeta*  en  again 
be  aeparated.  Bie  agnation  and  bv  )daxy  eondltlona  for  the  curvature 
tana  are 


%;ao)  =  0 

/\j  oo  ^ 

Tor  the  allp  and  teaverature-Juagp  effeeta 

\X  ^4o  X 

to)  = 


(7.22) 


1 


^  £>0 


0 


(7.23) 


fiw  ablution  of  (7.23)  can  be  written  down  bgr  Inapeetlon 


^  ^  K  (0)  0-  \) 


(7.24) 
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ItaA  eqaatlai  and  touadaxy  eoadltloaa  for  tha  dlaplae— nt  offoet 


term  are 


2. 


vhieh  bM  the  eolutloii 


I 

0 


Tioellyi  for  the  vortleity  ten 

^  ( ^>>/)  =  "  4"' 


(F.25) 


(F.26) 


'O'*  ‘  ^ 

7-*"  0 

Bie  entire  te^perature-funotlon  correction  term  can  then 
be  written 

^ 0-  ^  ^  *■  ^ v.p.v^,v| 

(F.aB) 

Sils  conpletes  the  preaentation  of  the  differential  eqaationa 
and  boundaxy  eondltions  for  the  conatant- fluid-proper^  oaae.  It  la  of 
intereat  to  obaerve  that  these  conatant- fluid  propertgr  aolutiOBB  are 
actually  identical  to  the  Tariahle-fluidppraperty  ahlutiatta  for  the 
apedal  caae  of  W  »  I  .  aUa  la  ohaerrahle  by  la«peetiag  (2.66), 
and  notlim  that  1(4 \  ia  tha  adlution  for  thia  i^paelal  oaae. 
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Uliag  thla  xtiult  in  th*  eq^ntioiif  of  Chapter  n  reduces  these  eqiuations 
to  those  of  the  present  eppenUx. 


AFfinSX  0 


lltetaixi*-Coapatation  Frooadiir* 

Hm  differential  eqioatloau  elth  'bouadexy  eoodltloiui  that 
were  solved  by  the  Burrouc^  220  eong^ter  eea  be  deeorlbed  as  a 
tw^poiat  bouadaxywvalue  problea.  For  the  bouadarT^layer  oaae 
the  eq^atlaas  are  nonlinear >  for  the  eorreetlon  texas  linear. 

Hw  eosouter's  flrat>Qarder  alaultanecnia- differential- eqyiatloas 
routine  (vhleb  uses  the  ftinge-Eutta  aethod)  was  applied  to  the 
problea.  Oiree  pairs  of  starting  values  (for  (o)  and  ^  (<>) 
wore  asauned,  and  then  the  (Invlseld  asyaptotle)  behavior  of  the 
functions  at  ^  —  7  was  used  to  obtain  successively  better  aad 
final  starting  values  by  means  of  double  lntezpolatlon>  whldi  was 
prograaned  on  the  aachlne.  Xhls  procedure  was  then  successively 
repeated  with  the  decreaents  of  halved,  until  the  final 

staurtlng  values  did  not  change. 
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APFflOn  1 

CcapariMB  of  Ihaoxy  vlth  Ii^ezlMBBts  of  Tmrfik  and 

Sie  tvo  roferanees  that  daserlte  ttao  above  ajqpaxlaaatB  are 
(to)  and  (4l);  la  this  appendix  they  vUl  be  referred  to  as  Barta  1 
and  II  reapeotlvely. .  Oae  foUcwlag  ealeulation  procedure  was  uaed 
to  coapare  the  present  theory  with  the  ejqperiaeatal  Husselt  nind)ers 
presented  In  Table  H,  p.  Bart  I.  Blrst>  the  velocity  gradleiit 
at  the  sta^EBatlon  point  had  to  be  detendned.  This  eas  done  by  using 
the  e^gperlBentally  detexaiaed  values  given  in  Table  a6,  Bart  II. 
Using  the  deflnltlooe  given  on  p.  9,  Bart  n,  the  relation  between 
the  aoaenclature  of  the  present  analysis,  and  of  Tewflk  and  QLedt 
is  the  fOUoMing 


AR 

V 


(H.1) 


where  €.  is  the  density  ratio  across  the  shock.  In  the  Nae)^ 

number  range  considered  (1.3)  to  3*7}>  and  around  free-Btrean 

pressures  of  10~^  ataosidiereB  (in  agreeaeent  with  the  "low”  stagnation 

(l8) 

pressures  of  80-120  microns,  asatloned  on  p.  14  nurt  I)  Kaufkan'  ' 
iadlcates  that  the  perfect-gas  formulae  are  sufficiently  accurate  to 
calculate  the  flow.  Bor  air  •■f  »  1'^  can  be  assumed,  and  then 
e  and  the  Mach  number  after  the  shock  can  be  obtained  from 


M—  Forthexmore  after  the  shock  an  isentropic  cos^ression 
takes  place,  for  which  ^  T  a  T  ,  and  the  variation 

at  viscosl-^  with  temperature  can  be  assumed  to  be  ^  > 

Then  the  Bqmdlds  nukber  that  is  of  interest  in  the  present  analysis 
can  be  related  to  given  by  Tavflk  and  dladt  as  follows: 
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-  -  A-SHteii  _ 

ii  -  ^2T?.V 


(H.2) 


Table  V  of  Bert  I  Showe  that  the  recovery  factor  le  almye 
elnost  exactly  1 .0  at  the  stagnation  point,  indicating  that  the 
adiabatic  nan  teaqierature,  T^n/  ,  ebich  foxaw  the  baale  of 
the  e^qieilaental  Bhaselt  xnabers  that  are  preaented,  ia  alaoat 
exactly  the  aaae  aa  the  ateguatlon  teagperature,  ,  etaldi 

fozaa  the  baala  of  the  heat-tranafer  fonulae  of  the  preaent 
analyaia.  Sien  the  not^dlaenalonal  heat-tmnafer  rate  of  the 
preaent  analyaia  can  be  related,  to  of  Tevfik  and  OLedt 

aa  foUom: 


‘  (H.3) 

etaere  it  eaa  aaauaed  that  during  the  laentropic  coigpreaalon  frcn  the 

I  -^o-s? 

ahock  to  the  stagnation  condition  \  Bie  eaqperlnental 

reault  obtained  in  (H.3)  can  then  be  coepared  to  the  theory  of  the 
preaent  analyaia,  aa  given  in  Chapter  lY.  The  dlq?laeeaent-eorreetlon 
ten  Hill  not  have  to  be  conaldered  beoauae  the  expexlnentally  aeanured 


velocity  gradlenta  already  Include  thla  effect.  The  paraaeter  defining 
the  nagnltude  of  the  curvature  correction  ia  the  inverae  aqioare  root 
of  the  Baynolda  number  given  in  (1.2).  FOr  the  allp  and  tenperatore- 
Jukp  corraotion  ten,  the  paxaaeter  that  ia  algnlfleant  la  the  ntlo 
of  the  nean  flee  path  at  the  vail  to  the  boundary-layer  thidkneaa*. 
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^  U  s»[«^  i  ra  ^U  K  2^Jt,j«7;i  ’ 


[T+oThJ" 


(lA) 


Teagpenture  ratio 


W 


In  thn  above  eagpresaiona  is 


tabulated  In  Table  A^,  p.  36,  Part  n.  Vlaeoally  ratio  ^ 


eaa  related  to  W 
.(10) 


M-rC 


by  using  the  viscosity  vs  teaqperature  varlatlan 
given  by  Haaaen^"''  and  (at  the  low  tesgperaturea)  by  Jakcb^^’^\  Por 
a  stagnation  tenperature  of  300 IC  (p  3,  Bart  l)  a  good 

approxlaatlon  for  this  vlaeosity  variation  is  /c  T° 
near  W'O-'b  near  W'0.7 

low  all  quantities  in  (eA)  are  known  except  K ,  .A  reasonable 
guesa  (e.g.  references  (36)  and  (27)  )  is 


k-|f 


Hi. 


(  v.t.  \  ) 


=  \-d 

'X. 


(H.5) 


How  all  quantities  tbat  are  necessary  to  apply  the 
theoretical  results  of  Qiapter  lY  hanre  been  deteained.  Using  the 
results  of  this  theory  (Table  Ill)>  the  hea'^transfer  rate  at  the 
stagnation  point  can  be  vriten  for  the  eonstant-fluiiUproperty  ease 

C\'W)  (}.».+ 


=  O.SUI o-'iS  - 


(H.6) 
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n*  eomipoBdlag  waibers  for  tint  Tulol^fluldpproportgr 
cue  eaa  be  detezalaed  as  fuactlonsj  of  W  uslag  the  gnqphe  of 
flgart  7.  Tbe  remits  of  tbese  ealeulatloBS  ere  presented  la  tezas 
of  tbe  difference  between  tbe  esletaleted  result  end  tbe  heet-trsasfer 
rate  based  on  eonstsat-fluld>»property  boundarjr-lsgrer  tbeoiy,  as  a 
fraction  of  tbe  latter.  The  experlnsntal  results^  as  eeloulated  la 
(b.3)i  are  presented  In  the  seas  aanner,  e.g. 

4  ~  Q,. 

(H.7) 

etc.  TlnaUy  dlsplaeeaent  effect  paxwaeter  Is  also  tabulated 

for  tbe  sake  of  reference.  Tbe  calculation  of  It  Is  based  on  an 
snplrleal  relatlan  between  stagaatlonppaiat  velocity  gradients 
for  circular  cylinders  and  Naeb  nuBiber>  as  presented  by  Besbotko  and 
Beckwltb^^^^ .  (This  is  tbe  ''lnfinlte''-BByaolds>nuBd)er  ease).  UOlag 
these  and  tbe  experlswntally  detezalaad  velocity  gradients  of  TOwflk 
and  (Uedt  (H.l)  la  e^^resslon  (2.22)  Is  detendjsed  as 

follows 


^5 


A 


A 


-  A 


Wait  wV tv 


(H.8) 


IGNDraATOIlE 


a.  speed  of  souBd 

(with  mbeerlpt)  constant  coefficients  (Cbapter  I  only) 
/V  stagnation-point  velocity  gradient;  u«Ax 

(r  constant  coefficients  (Oiapter  1) 

Cj  displacenent-effect  paranster  (equation  2.22) 

(-S  coBdiiaations  of  Bessel  functloou  (equation  3*^) 

Cf  specific  heat  at  constant  pressure 


e  density  ratio  across  (nomal)  shock; 

a,  vectors  in  orthogonal  systssi  (Appendix  A) 

Ec  )  energy- equation  differential  operator  (equation  2. to) 

)  inhoangeneous  terns  in  energy  equation  (eqjuation  2. to) 

functions  in  esqpansion  of  viscous  flow  about  stagnation  pointy 
X*  and  xt  tezns  respectively.  'Without  adscript  stream 
function,  vlth  adscripts  pressure,  density, 

and  temperature  respectively.  (S.g.  is  teaiperature- function, 

term,  etc.,  equations  2.2  through  2.$) 

V\  enthalpy 

metric  functions  in  orthogonal  coordinate  system  (Appendix  A) 

H  altitude,  ft. 

X,C  )  Bessel  function  of  C  )  (Chapter  Ill  and  Appendix  C) 

4.  heat  conductivity 

constant  of  integration  (eqviatlon  2.33) 

K,  ,  Kf.  proportionality  constants  (equation  2.21) 

(  )  Bessel  function  of  (  ^  (Chapter  IH  and  Appendix  C) 

L  reference  length  indicative  of  body  size,  detexmlnss  ^  ^  °  X 
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H  Nach  nmAMr 

M(  1  maMntuifr>«gu«tioii  dlfferwitlal  operator  (eqiuatlon  2.1(0) 
M(.C  ^  lahonogeaeous  terms  In  momentum  equation  (equation  2.k0) 
Y\  0  tor  tvD-diaensicnal,  1  for  axially  eyometrle  flov 

Huv/  Niisselt  number 

Cr(  Order  of  (.  ) 

p  presstire 

Prandtl  number 


Q 

r 

K 

d 

0,vr 

V 

V 
w 

^  ,  '5r 

o(. 


velocity  vector 

heat- transfer  rate  normal  to  surface  ^ 

radial  coordinate 

nose  radius  of  curvature 

gas  constant 

Beynolds  number 

absolute  temperature 

velocity  ccnponents  in  <  and  directions  respectively 

free- stream  veloci'ty;  ft/sec 

Vortlclty  parameter  for  axially  symnetzle  stagnatlonrpolnt 
flow  (equation  1.19) 

cooling  ratio;  W 

boundary-layer  coordinates  (equation  1.1) 

^  mC 

exponent  of  temperature  for  specific  heat;  \ 

aecoBudation  coefficient  for  energy  transfer  at  solldpgas 
interface 

ratio  of  specific  heats 

boundary- liyer  thickness  parameter  X  -  |^ 


r 

i 

a 


Inviscld  shock  standpoff  distance 
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£  •xpcBUt  Of  taaporature  foor  boat  eondnetlTltgr;  ^  > 

ft!  bOKBAuy^Iigror  eoordlait*  (•foatloa  2.1)  />?  • 

diaplaecaent  maiber  (•qqatloii  2.3if) 

^  (vttboat  auiboexlpt)  eoovdliiata 

^  (idtb  •ubsczlpts)  taaperature  fuaetion  la  eoBStaat^flaldp 
property  flow  (eguatlone  7.7  and  F.20) 

A  aeea  free  path  la  gae 

vlseoclty 

y  klaeaatle  vieeoelty 
f>  densll^ 

^  fraetloo  of  dlffoeely  reflected  aoleeulee 
T  abear  parallel  to  surface  X*  4^. 

coapressible  stream  fuactlom 
^  dlssipatloa  fooetlom  la  energy  egoatloa 
ou  exponeot  of  temperature  for  viscosl'^  ^  «<  I 
SI  Torticlty  SI  » 
proportional  to 

Subaezipts; 

^  p.„suecesslTe  terns  la  nyaolde-ataiber  expansion;  bouadazy^lagrer 
^  1st  order,  2ad  order,  etc.  coarzection  terms;  also  (la  CSmipter  I) 
successlTe  coefflciaats  la  u  expaasioa  of  lavlsoid 

quaatltlee.  ^ 

c  eozTeetlon  due  to  curvature  effect 

com^  result  of  eonpresslUe  (Tariable-fluld>property)  analysis 
oorrectlQn  due  to  dlsplaeonsat  effect 
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t/xp  expariJMntAl  remit 

i  correction  due  to  teapemture  Jiaq^  C  K ,  ^  K&) 

Cw('  result  of  ineoaq^reseible  (constent-fluidpproperty)  enelysls 

M  correction  due  to  finite  asea  free  peth  (eqpiatlan  2.39) 

P  refers  to  pressure 

r  refers  to  density 

r  pertlal  derivative  vlth  respect  to  r  (Appendix  A  only) 

reference  quantities  (see  Chapter  IV)>  equations  (4.1)  end  (4.3) 
(except  (F.23)  sad  (7.24). 
s  Invlscld  sta^atlon  value 

condition  after  shock 

SL  correction  due  to  velocity  slip 

st-’.i/ftL  standard  conditions  at  standard  ataosphere 

t  refers  to  tenqwrature 

V  correction  due  to  vortlclty  effect 

W  condition  at  solid  surface  ("wall") 

partial  derivatives  with  respect  to  x  and 

^  partial  derivative  with  respect  to  ^ 

free- stream 

Siverscrlpts 

'  derivative  with  respect  to  ^ 

» 

derivatives  with  respect  to  temperature 
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ADDENDUM 
February  I962 

bir.ae  completion  of  the  research  work  presented  in  the  foregoing 

report  several  new  developments  in  low- Reynolds -number  boundary-layer 

(A1,A2) 

theory  came  to  the  author's  attention.  In  particular  Vandyke  and 

(A3) 

.^aslen  utiliaei  a  Lagerstron-Cole  type  expansion  procedure  to  calculate 

eeconi-orier  terms  (i.e.  first-order  correction  terms)  to  stagnation-point 

and  other  boundary- layer  problems.  The  purpose  of  this  Addendum  is  to 

clarify  a  shortcoming  of  the  foregoing  analysis  which  became  apparent  in 

cot.parison  with  these  new  developments,  ana  which  was  brought  to  the 

author's  attention  by  Professor  Vandyke  in  private  conversation.  In 

addition,  an  important  detail  regarding  the  application  of  the  foregoing 

theory  to  experimental  low-iteynolds-number  stagnation-point  measurements 

will  also  be  clarified;  this  was  brought  to  the  .  'thor's  attention  by 
(Aa) 

Prof.  Rott  in  a  differ-nt  context. 

net  the  external  flovf  be  described  by  a  streamfunction  of  the 
form;  ‘ ~ 


where  the  equation  above  replaces  (1.22).  The  first-order  correction 
terms  appearing  in  the  above  expression  are  not  necessarily  related  to 
the  nose  radius,  rather  (as  explained  in  the  introduction  to  Chapter  II.) 
x  is  merely  a  convenient  reference  length.  The  pressure,  temperature, 
and  density  expansions  can  be  related  to  the  streamfunction  by  means  of 
the  procedure  described  in  Chapter  1.  (This  is  permissible  for  the  first- 
order  correction  tenriS  because  the  viscous  terms  in  the  equaticxis  axe  of 
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February  1962 

6ir.ce  completion  of  the  research  work  presented  3n  the  foregoing 

report  several  new  developments  in  Icw-HeynolJs-nnmber  boundary-layer 

(A1,A2) 

theory  came  to  the  author's  attention,  in  particular  VanJyke  and 

(A3) 

iaslen  utilized  a  Lagerstrom-Cole  type  expansion  procedure  to  calculate 

Eeconi-orJer  terms  (i.e.  first-order  correction  terms)  to  stagnation-point 

and  other  boundary- layer  problems.  The  purpose  of  this  Addendum  is  to 

clarify  a  shortcoming  of  the  foregoing  analysis  which  became  apparent  in 

coE.parison  with  these  new  developments,  aiu  which  was  brought  to  the 

author's  attention  by  Professor  ianJyke  in  private  conversation.  In 

aidition,  an  important  detail  regarding  the  application  of  the  foregoing 

theory  to  experimental  low-iteynolds-number  stagnation-point  measurements 

will  also  be  clarified;  this  was  brought  to  the  author's  attention  by 
(Au) 

Prof.  Hott  in  a  different  context. 

Let  the  external  flov/  be  described  by  a  streamf unction  of  the 


(Ad.l) 


where  the  equation  above  replaces  (1.22).  The  first-order  correction 
terms  appearing  in  the  above  expression  are  not  necessarily  related  to 
the  nose  raiius,  rather  (as  explained  in  the  introduction  to  Chapter  II.) 
u  is  merely  a  convenient  reference  length.  The  pressure,  temperature, 
and  density  expansions  can  be  related  to  the  streamfunction  by  means  of 
the  procedure  described  in  Chapter  1.  (This  is  permissible  for  the  first- 
order  correction  tenris  because  the  viscous  terms  in  the  equations  are  of 


order  —  . ) 


Only  the  pressure  expansion  is  of  interest  here; 


The  "inner"  expansions  of  the  saire  variables  are  as  given  in 
Chapter  11.  litching  the  two  exjiansions  requires  that  as  the 

"inner"  expansion  approach  the  external  flow  described  by  (Ad.l)  and 
(Ad. 2).  This  implies  that  the  solution  for  f;^  ,  f  ,  f  ,  and 

fp^  are  as  given  in  Chapter  li.  Furthermore,  usirig  the  boundary-layer 
solution  in  the  Hatching  requirement  one  obtains  that; 


(Ad. 3) 


ainoe  any  constant  that  might  appear  in  the  limiting  behaviour 
of  (i^)  or  higher-order  corrections  would  be  of  order  ^  or  smaller, 
it  is  apparent  that  the  constant  a„  in  (Ai.l)  can  be  identified  by 
comparison  with  (Ad. 3)  as; 

■L.  =  -  K.y* 

It  is  likewise  apparent  that  to  determine  constant  bg  the 
limiting  behaviour  of  f,  (^)  must  be  known.  However  equaticxi  (Ad. 2) 
shows  that  b||  appears  as  a  boundary  condition  necessary  for  the 
integration  of  gp,  ((iy),  which  in  turn  is  necessary  for  the  solutitai  of 
f|(y).  I^is  shows  that  b,|  can  not  be  determined  within  the  framework 
of  stagnation-point  flow  alone;  it  can  then  be  defined  as  the  undeter¬ 
mined  "lisplacement  constant",  i.e.; 


(For  further  discussion  on  the  somewhat  aiibiguous  nature  of 


defining  the  displacement  effect  reference  is  made  to  (Au).) 

Using  results  (Ad.U)  and  (Ad.^)  in  (Ad.l)  and  (Ad. 2)  the 
'•outer** -flow  expressions  become; 


H'*  ...  -Kijl  (3>'jr +•']>•• 


(Ad.  6) 


f  =  Pi-s/ 


(1+nf-  I  d-  ^  +...  + 


Corresponding  to  (Ad. 6)  the  constant  of  integration  in 
integrating  gp|  (a^)  (of.  pp.  31“32)  becomes; 


(Ad. 7) 


which  now  replaces  (2.35).  Then  the  corrected  form  of  (2.36)  becomes; 

0+"vf+'^ \i,*-CV^- Q>-  (Ad, 


s,A- 

Thus  it  is  clear  that  the  momentum  equation  for  the  vorticity 
correction  term  is  not  hranogeneous  as  given  in  (2.57)»  but  rather; 


r« 


(Ad. 9) 


is  the  correct  equation.  Now  the  behaviour  of  the  momentum  equation  for 
the  correction  tera  as  can  be  observed; 


iliuB  it  is  clear  from  (2.57)  that} 

Un  the  other  hand  (Ai.6)  and  (Ai.?)  imply  that; 
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(Ad. 11) 


(Ad.  12) 


Finally,  the  limiting  behaviour  of  the  displacement  correc¬ 
tion  term  consistent  with  (2.55)>  (Ad. 6)  and  (Ad. 10)  is; 

^  (2.53) 


observing  that  the  energy  equation  and  all  the  remaining 
oo\)niary  conditions  for  the  displacement  and  vorticity  correction  terms 
are  identical; 


—  0  (Aa.l3) 


and  using  the  five  foregoing  expressions,  it  is  easy  to  see  that; 

satisfies  both  the  correct  momentum  equation  (Ad.?)  and  boundary 
condition  (Ad. 12).  The  energy  equation  and  the  remaining  boundary  con¬ 
ditions  are  automatically  satisfied  by  (Ad.lu).  Equation  (Ad.lU)  is 
then  the  correct  vorticity  term;  the  temperature  function  ft,y^^^  is 
correspondingly  modified.  The  correction  to  the  vorticity  effect 

described  by  (Ad.lU)  has  become  known  in  the  literature  as  the  "vorticity- 

(A5)  (Al) 

inauced  pressure-gradient"  effect  (e.g.  Li  ,  Vandyke  ,  etc.); 
and  the  present  result  and  the  results  of  the  above  mentioned  investi¬ 
gations  are  thereby  in  agreement. 
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One  may  now  consider  an  application  of  the  foregoing  theoiy 
to  experimental  measurements  of  low-Reynolda-riumber  stagnation-point 
flow.  Let  the  velocity  gradient  A  be  based  on  pressure  distribution 
data  obtained  at  a  "low"  Reynolds  number,  so  thatj 


a'  - 


^•0 


(Ad. 15) 


Using  (Ad. 8)  and  (2.5)  in  the  above  expression  one  obtains; 


(Ad. 16) 


i/-,  “  -  ♦  I  f.  =  |l  - 

where  A  corresponds  to  an  "infinite"  Reynolds  number. 

»Now,  considering  only  the  highest-order  corrections,  A  can  be 

expressed  in  terms  of  A  ; 

exp 

where  constant  K^,  is  defined  by; 


(Aa.l?) 


(Ad. 18) 


if  now  A  is  to  be  used  to  predict  physical  quantities  in 
exp 

the  boundary  layer,  and  if  the  Reynolds  number  is  sufficiently  low  to 
necessitate  consideration  of  the  first-order  correction  terms,  then  the 
change  described  by  (Ad. 17)  and  the  effect  of  this  change  on  the  boundary-1  eyar 
term  must  also  be  included  for  a  consistent  theoretical  prediction  to  this 
order.  The  effect  on  the  boundazy-layer  quantities  of  changing  A  to  the 
first  order  is  already  known  (of.  pp.  38-39);  it  is  the  displacement- 
effect  term.  Thus  expressions  of  the  form  given  in  Chapter  iV., 
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e.g.  (ii.6)j 

.u  ,  y*  ♦  i  ji 


•1-  nV 


must  now  be  revise:,  in  accordance  with  (Ad.l?); 


H  s  y*  4  -  ] 

A 

Uhv}^ 

^  Wlw 

1  '1 

.  A  f  / 

U.3,  ^ 

A  f 

^  fVA  L 

fv» 

i  A. 

-  y*  (Ai.20) 

VA  U^^i. 


The  two  expressions  above  differ  because  (Ad.l?)  uses  a  velocity 
graaient  A  calculated  from  inviscid  theory,  or  reasured  at  an  equivalent 
"infinite"  (i.e.  very  high)  Heynolds  number,  whereas  (Ad.20)  is  based  on 
A^^p  obtained  from  pressure  measurements  at  the  particular  ("low") 

Reynolds  number  corresponding  to  .  The  former  method  must  therefore 
account  for  an  undetermined  change  in  A  due  to  uhe  displacement  effect, 
whereas  the  latter  has  already  taken  this  into  account,  but  must  also 
account  for  the  changes  in  surface  pressure  gradient  due  to  the  centrif¬ 
ugal  pressure  rise  and  vorticity  interaction  effects. 

Constant  (Ad.l8)  is  a  function  of  M  and  is  tabulated  below, 
as  a  supplement  to  Table  111. 
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TABLE  Ad.  I.  Values  of  K^,. 


H 

K 

C 

n=C 

n»l 

1.0 

“DTW32 

u.blbb 

0.75 

0,7998 

U.6927 

0.5 

u.  dliOit 

0.5601 

0.25 

0.a820 

0.dl60 

0.1 

0.3717 

0.3216 

The  experimental  comparison  presented  in  Chapter  IV.  and  Table 
iV.  and  carried  out  in  Appendix  H  must  also  be  revised  in  accordance  with 
(Ai.20).  Values  of  and  •ijjCO)  clearly  indicate  that  the  theoretical 
heat-transfer  rate  will  be  somewhat  higher  than  the  predictions  pre¬ 
sented  in  Table  IV,,  thereby  making  the  agreement  between  experiment 
and  theory  less  favorable  than  implied  by  that  Table. 
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Subscripts : 

Gorr  corrected  vorticity  tern  due  to  induced  pressure  effect 
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